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SECTION  I 


INTRODUCTION 


The  research  effort,  conducted  during  the  last  years  at  the  Department  of 
Aeronautical  Engineering  at  the  Technion,  developing  improved  models 
for  air  combat  analysis,  has  shown  that  an  analytical  investiga¬ 
tion  of  a  complex  dynamic  problem  can  lead  to  meaningful  and  applicable  results. 
USAF  support  of  the  research  program  guaranteed  that  the  effort  be  focused  on 
relevant  air  combat  problems  and  made  possible  to  complete  the  investigation  in 
a  relatively  short  time. 

It  has  been  a  concensus  for  many  years  that  the  conflicting  nature  of  air 
combat  should  be  analyzed  by  a  differential  game  approach.  However,  the 
inherent  mathematical  complexity  of  such  game  formulation  with  realistic  dynamic 
models  left  little  hope  for  a  comprehensive  investigation.  In  order  to  be  able 
to  obtain  practically  useful  results  it  was  first  necessary  to  identify  an 
analytical  method  of  approximation  which  enables  some  order  reduction  of  the 
mathematical  model.  The  technique  of  multiple  time  scale  singular  perturbation 
seemed  to  be  an  attractive  candidate  based  on  its  successful  application  in 
aircraft  performance  optimization  problems. 


The  first  phase  of  the  research  program  was  oriented  to  overcome  the 
conceptual  difficulty  in  applying  the  singular  perturbation  approximation  to 
nonlinear  zero-sum  differential  games  [1],  In  the  sequel,  a  set  of  examples  of 
gradually  increasing  complexity  were  solved  leading  to  development  oi  a  consistent 
methodology.  This  methodology  provides  a  uniformly  valid  zero-order  solution  in 
a  feedback  form  for  the  control  strategies  of  a  realistic  three-dimensional 
variable  speed  pursuit -evasion  game.  A  basic  operational  scenario  of  future 
air  warfare,  the  medium  range  air-to-air  interception  using  advanced  homing 
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missiles,  serves  as  a  suitable  example  to  demonstrate  the  merits  of  this 
innovative  approach. 

The  major  effort  of  the  investigation  in  the  last  year  (sponsored  by 
AFWAL  Contract  No.  F33615-81-K-3007)  was  concentrated  on  developping  a 
computer  program.  This  program  implements  the  variable  modelling  forced  singular 
perturbation  technique  to  solve  the  problem  of  medium  range  interception  (using 
missiles)  formulated  as  a  perfect  information  zero-sum  differential  game. 

The  objective  of  this  report  is  to  document  the  results  obtained  during 
this  last  year  and  to  outline  the  potential  of  the  methodology  developed  for 
future  applications.  The  report  starts  (Ch,  2)  by  summarizing  the  results 
obtained  in  the  previous  phases  of  the  multi-year  research  effort,  followed 
(Ch.  3)  by  the  formulation  of  the  three-dimensional  medium  range  time  optimal 
air-to-air  interception  game.  The  solution  of  this  problem,  presented  in  Ch.  6, 
is  reached  as  a  consequence  of  modelling  considerations  (described  in  Ch.  4) 
and  the  insight  gained  by  solving  first  (in  Ch.  5)  the  interception  game  in  the 
vertical  plane.  In  Ch.  7  an  imaginary  air  defence  scenario  serves  as  an  example 
to  demonstrate  the  usefulness  of  the  computationally  efficient  solution  technique. 
In  the  following  chapters  the  accuracy  of  the  zero-order  approximation  is 
analysed  and  eventual  corrections  are  discussed.  The  potential  merits  of 
applying  the  zero-order  feedback  solution  for  airborne  implementation,  as  well 
as  for  rapid  systematic  performance  analysis,  are  presented  in  the  concluding 
chapter. 


2 


SECTION  II 


SUMMARY  OF  PRECEDING  STUDIES 

1.  APPLICATION  OF  SPT  FOR  NONLINEAR  DIFFERENTIAL  GAMES 

The  application  of  singular  perturbation  technique  (SPT)  for  solving 
certain  types  of  pursuit -evasion  problems  formulated  as  non-linear,  zero- 
sum  differential  games  has  been  demonstrated  in  the  past  [1].  Motivation 
for  using  this  approximation  technique  was  explained  by  the  necessity  to 
overcome  the  high  dimensionality  (as  well  as  the  non-linear  character)  of 
realistic  pursuit -evasion  problems. 

SPT  is  based  on  the  assumption  that  some  components  of  the  state  vector 
behave  as  "fast",  compared  to  the  other  "slow"  ones,  leading  to  an  eventual 
time  scale  separation.  Using  this  approach,  the  original  problem  is 
separated  into  several  subgames  of  lower  dimensions,  which  can  be  solved 
analytically.  This  process  is  meaningful  only  if  it  is  assumed  that  the 
original  game  has  indeed  a  saddle  point  solution  (an  assumption  which  un¬ 
fortunately  cannot  be  verified). 

Dealing  with  pursuit-evasion  games,  the  reduced  order  game  represents  a 
"simpxe"  pursuit,  where  the  "fast"  state  variables,  i.e.,  flight  direction  and 
velocities  of  both  opponents,  are  considered  as  controls.  The  reduced  order 
game  has  a  smooth  value  function  satisfying  the  sufficiency  condition  for 
saddle  point  optimality  [1] . 

This  game  is  unable  to  satisfy  the  initial  and/or  terminal  conditions  im¬ 
posed  in  the  original  problem  on  those  variables  for  which  the  dynamics  is 
neglected  in  the  reduced  game. 

In  order  to  match  these  conditions  several  boundary  layer  games  are 
required.  The  boundary  layer  games  are  usually  decoupled  into  one-sided 
optimal  control  problems,  for  which  satisfaction  of  "game  optimality"  is  not 
required.  Based  on  these  subgames  a  uniformly  valid  composite  solution  is 
synthesized  and  it  represents  a  suboptimal  zero-order  approximation  for  the 
exact  solution  of  the  original  game  (if  such  solution  indeed  exists) . 


The  outcome  of  the  original  game  using  the  suboptimal  SPT  strategy  pair  is 
defined  as:  Extended  Value  of  the  game.  It  was  shown  [1]  that  this  extended 
value  is  bounded  on  both  sides  satisfying  a  week  saddle  inequality.  Further¬ 
more,  it  was  proven  [1]  that  if  the  perturbation  parameter  (e)  becomes  very 
small  the  Extended  Value  of  a  singularly  perturbed,  zero-sum  differential  game 
approaches  as  a  limit  the  exact  Value  of  the  game. 


2.  SPT  SOLUTION  OF  THE  "GAME  OF  TWO  CARS" 

The  usefulness  of  SPT  has  been  demonstrated,  first,  via  the  example  of  the 
"Game  of  Two  Cars"  (a  2-D,  constant  speed  pursuit-evasion  game) .  In  this 
example,  the  original  equations  expressed  in  polar  coordinates,  were  transformed 
in  order  to  discover  the  existence  of  a  small  parameter  (the  ratio  of  aircraft's 
turning  radius  to  the  initial  distance  of  separation  R/rg) •  As  a  result,  the 
singularly  perturbed  differential,  zero-sum  differential  game  had  a  simple 
zero-order  closed  form  solution.  The  strategy  of  each  player  is  composed  of 
three  phases: 

(a)  a  "hard"  initial  turn  until  the  velocity  vector  is  aligned  with  the  line 
of  sight; 

(b)  a  subarc  of  line  of  sight  guidance  if  the  other  player  is  still  in  phase  (a); 

j 

(c)  a  straight  line  dash  till  capture. 

This  approximation  was  compared  to  the  exact  solution  of  the  same  problem 
given  by  Simakova  [2], 

In  the  exact  solution  the  direction  of  the  final  dash  is  the  common  tangent 
of  the  turning  circles  of  the  players  as  determined  by  the  initial  conditions  of 
the  game.  The  optimal  strategy  of  the  players  is  to  align  their  velocity  vector 
with  this  tangent.  If  the  initial  conditions  are  such  that  both  players  finish 
their  hard  turn  simultaneously  the  trajectories  of  the  optimal  and  the  SPT 
solution,  and  consequently  the  outcome  of  the  game  are  identical.  For  any  other 
initial  conditions  one  of  the  players  completes  its  "hard"  turn  earlier  and 
following  the  suboptimal  strategy  of  phase  (b)  deviates  from  the  optimal  one. 
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If  it  is  the  pursuer,  the  time  of  capture  will  be  slightly  longer  than  in  the 
optimal  game.  If  the  evader  reaches  first  the  line  of  sight  direction,  the 
capture  time  of  the  SPT  solution  will  be  shorter  than  the  Value  predicted  by 
Simakova.  The  differences,  however,  are  small. 

A  quantitative  example  shows  a  remarkably  good  agreement  between  the  outcome 
of  this  game  using  suboptimal  strategies  and  the  exact  Value.  Even  for 
relatively  large  values  of  the  perturbation  parameter  (e  =  0.4)  the  difference 
is  only  a  few  percent  (the  relative  error  in  this  example  is  between  the  order 
of  e2  and  e3) . 

3.  VARIABLE  SPEED  MEDIUM  RANGE  INTERCEPTION  GAME 
3.1  Model 

As  a  second  step,  a  variable  speed  horizontal  pursuit-evasion  game  of  two 
airplanes,  with  realistic  thrust  and  aerodynamic  data  has  been  analysed  [3]. 

In  this  case  it  was  difficult  to  show  the  existence  of  a  small  parameter  of 
physical  significance  which  can  transform  the  system  into  a  singularly  perturbed 
structure.  However,  time  scale  separation  of  the  variables  was  known  to  exist, 
therefore,  a  small  parameter  inserted  artificially  into  the  system,  multiplying 
the  derivatives  of  the  "fast"  variables. 

This  so-called  Forced  Singular  Perturbation  Technique  (FSPT) ,  has  been  used 
extensively  for  aircraft  performance  optimization  [4-7],  FSPT  was  formally 
justified  [8]  as  a  legitimate  technique  to  produce  zero-order  feedback  solution 
for  problems  which  exhibit  genuine  time  scale  separation. 

Results  of  the  constant  speed  model  have  shown,  that  if  the  initial  separation 
is  large  compared  to  the  turning  radius  of  the  players  the  line  of  sight  dynamics 
is  very  "slow".  Eventual  speed  variations  are  assumed  to  be  next  in  the 
hierarchy  while  the  directions  of  the  players  are  considered  as  the  "fastest" 
variables . 

Assuming  that  capture  range  is  greater  than  the  pursuer's  turning  radius,  a 
uniformly  valid  zero-order  approximation  of  the  optimal  strategies  has  been 
obtained  analytically  in  a  feedback  form. 


The  solution  indicated  that  for  large  initial  separation  the  optimal 
strategy  terminates  by  a  "tail  chase"  flown  at  maximum  velocities.  This 
reduced  order  solution  is  reached  by  both  airplanes  by  performing  a  variable 
speed,  full  thrust,  turning  maneuver  determined  by  a  time  optimal  compromise 
between  the  lateral  and  longitudinal  accelerations.  Implementation  of  the 
suboptimal  feedback  strategies  requires  to  measure  only  the  relative  angular 
position  and  the  vehicle's  own  velocity.  It  seems,  therefore,  very  attractive 
lor  real  time  on-board  applications. 

Numerical  example  solved  in  [3]  confirmed  the  assumptions  of  time  scale 
separation  justifying  the  forced  singular  perturbation  approach. 

In  this  classical  FSPT  approach  the  hierarchy  of  the  state  variables  (their 
respective  role  as  being  "slow"  or  "fast")  is  preserved  for  the  whole  duration 
of  the  game.  Experience  has  shown,  however,  that  in  some  nonlinear  dynamic 
problems  the  relative  rate  of  change  of  the  variables  may  vary  considerably 
in  the  course  of  the  solution. 

3.2  Model  B. 

The  inherent  property  of  the  previous  model  (defined  as  model  A)  was  the 
requirement  for  full  thrust  as  a  control  strategy  in  all  situations.  Engineering 
insight  indicates  that  such  a  model  may  not  be  representative  in  a  case  where 
one  of  the  airplanes  flies  at  a  speed  largely  exceeding  its  own  "corner  velocity" 
(where  the  instantaneous  turning  rate  is  maximum)  if  its  flight  direction 
deviates  considerably  from  the  one  required  for  optimality.  It  seems  intuitively 
clear  that  such  an  airplane  has  to  reduce  its  speed  (to  improve  its  turning  rate) 
as  rapidly  as  possible  which  implies,  for  horizontal  flight,  zero  thrust  as 
optimal  control  strategy.  In  this  case  speed  variations  cannot  be  more  considered 
being  "slow"  as  compared  to  the  constrained  turning  rate.  To  properly  model  such 
situation  it  has  been  proposed  [9]  to  declare  velocity  variations  as  the  "fastest" 
dynamics  in  the  SPT  hierarchy.  Turning  dynamics  were  still  considered  "faster" 
than  the  variations  of  the  line  of  sight  and  separation  distance.  This  new  SPT 
model  (defined  as  model  B) ,  allows  for  variable  throttle  setting  as  optimal 
control  strategy.  The  validity  of  this  model  is,  however,  conditional.  When  the 


model  dictates  full  thrust,  deceleration  rate  becomes  slower  than  turning  rate 
and  model  A  may  be  more  representative. 


-3.3  Combined  Strategy  Logic 

In  [9]  a  variable  modeling  approach  has  been  proposed.  This  model  is 
based  on  a  Combined  Strategy  Logic  (C.S.L.)  synthesized  by  using  both  A  and  B 
models.  The  first  step  of  C.S.L.  is  to  verify,  using  model  B  (for  each  player 
settarately)  whether  the  initial  conditions  of  the  engagement  require  a  zero 
thrust  subarc. 

A  negative  answer  leads  to  the  use  of  model  A  for  the  respective  player 
during  the  whole  game.  For  a  positive  answer,  model  B  is  used  until  full  thrust 
is  required.  Then  the  solution  continues  with  model  A.  Since  both  models  use 
the  same  reduced  order  game,  the  variable  modeling  computational  scheme  is  only 
slightly  more  complex  than  the  one  using  a  single  model. 


4.  THREE  DIMENSIONAL  CONSTANT  SPEED  MODEL 

Real  interception  is  a  three-dimensional  problem.  Therefore,  further 
extensions  of  the  previous  models  are  necessary.  As  a  first  step  towards  this 
goal,  the  2-D  constant  speed  model  (the  game  of  two  cars)  has  been  extended  to 
include  the  effects  of  a  3-D  flight  [10). 

Similarly  to  the  2-1),  constant  speed  model  it  was  assumed  that  turning  rates 
of  both  aircraft  are  much  faster  than  the  rate  of  change  of  the  line  of  sight. 
With  this  assumption  a  general  solution  (using  vectorial  notation)  has  been 
obtained.  The  solution  indicates  that  the  optimal  strategy  for  each  aircraft 
is  to  orient  its  angular  velocity  vector  in  a  direction  perpendicular  to  the 
plane  defined  by  the  line  of  sight  and  its  current  velocity  vector. 

The  results  of  this  strategy  were  the  following: 

(a)  During  the  first  stage  of  the  engagement  each  aircraft  moves  in  its  own 
plane  (defined  by  the  line  of  sight  and  its  velocity  vector). 
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(b)  As  soon  as  one  of  the  aircraft  coincide  its  velocity  vector  with  the  line 
of  sight,  it  keeps  moving  along  this  line. 

(c)  During  the  last  stage  both  aircraft  move  along  the  line  of  sight. 

One  of  the  consequences  of  this  solution  is  that  a  game  which  has 
progressed  into  a  planar  configuration  will  remain  in  it  till  termination.  A 
fact  supported  by  other  references  [11-13]. 


5.  A  CONCLUDING  REMARK 

The  results  summarized  in  this  chapter  have  served  as  guidelines  for  the 
formulation  and  an  appropriate  solution  of  the  real  air  combat  problem  of 
interest:  "the  3-D  medium  range  interception  game",  using  a  3-D  variable 
speed  FSPT  model. 
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SECTION  III 


FORMULATION  OF  THE  3-D  MEDIUM  RANGE  TIME  OPTIMAL, 

AIR-TO-AIR  INTERCEPTION  GAME 

1.  APPLICATION  OF  SPT  FOR  NONLINEAR  DIFFERENTIAL  GAMES 

The  motion  of  a  point  mass  lifting  vehicle  over  a  flat  non-rotating  earth 

assuming  symmetrical  flight,  is  governed  by  the  following  set  of  nonlinear 


ordinary  differential  equations: 

• 

X 

=  V  cos  y  cos  X 

(3.1) 

• 

y 

=  V  cos  y  sin  y 

(3.2) 

• 

h 

=  V  sin  Y 

(3.3) 

V  -  g  \ 

rC  Tmax('h,V)C0s(a  +  V  '  D  1 

(3.4) 

[  W  Sin  YJ 

rC  T  (h,V)  sin(a  +  B~)  +  L  ^ 

(3.5) 

^  ^  1  M  -  cub  yj 

v  -  £  i 

rf,  T  (h,V)  sin(a+PT) +L  1 

1  max  T  -in  u  ? 

(3.6) 

x  ‘  v  1 

[_  W  cos  Y  J 

m  =  -  j  Cs(h,V,T) 

(3.7) 

L  =  4  P(h)V2  S  C. (a.M) 

t.  L- 

(3.8) 

D  =  -  p(h)V2  S  Cd(M,Cl) 

(3.9) 

For  a  parabolic  drag  polar: 

Cd(M,Cl)  =  CD  (M)  +  K(M)  C2(a,M)  (3.10) 


and  aerodynamic  load  factor  defined  by: 


n 


(3.11) 


L  A  7  p(h)  ^  CL^a'^ 

=  w  w7s 

The  drag  force  can  be  expressed  as: 

D  =  DQ  +  n2  Dj  =  D(h,V,n)  (3.12) 

where  is  the  zero  lift  drag  and  Dj  the  induced  drag  in  level  flight: 


D 


I 


2 

KW 

j  p(h)  V2s 


(3.13) 


2.  CONTROL  CONSTRAINTS  . 

The  aircraft  controls  are  throttle  (£) ,  angle  of  attack  (a),  bank  angle  (y) 
Hn<*  thrust  tilt  angle  with  respect  to  the  aircraft  longitudinal  axis  (8^) . 
a  can  be  replaced  as  control  variable  by  n  (or  CL)  according  to  the 
following  rule: 

Below  "corner  velocity"  defined  as: 


2n  W 
max 

'C  “  l.p(h)  SC  (M^IJ 

max 


=  f- 

lp( 


(3.14) 


The  load  factor  is  constrained  by  the  maximum  lift  coefficient: 


n  £  n  (h,M)  =  j  p(h)  V2  S  CL  CM) /W 

max 


(3.15) 


Above  the  "corner  velocity"  the  load  factor  is  constrained  by  the  aircraft 
structural  properties : 


n  *  n 


max 


(3.16) 


Throttle  constraints  are: 


10 


* 


0  $  £  £  1  (3.17) 

which  are  equivalent  to  the  following  constraints  on  the  aircraft  thrust: 

°  *  T  «  T  (h,V)  (3.18) 

max 

P  .nk  ingle  is  constrained  by: 

-it  £  v  £  n  (3.19) 


The  thrust  tilt  angle  constraints  are: 

6t  *  B,p  4 

min  max 


(3.20) 


•3 .  STATE  CONSTRAINTS 

Minimum  altitude  limit: 

h  >  0 


Maximum  dynamic  pressure  limit: 


1 ' 1  »<h>v2  ‘  Vx 


Maximum  Mach  number  limit: 


V  <  a(h)  M 


max 


Loft  ceiling  (h^)  limit 


1  ,,  ,  ,,2  .  W/S 

2  p(hL}  V  CT  (mT 


max 


(3.21) 


(3.22) 


(3.23) 


(3.24) 
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4.  SIMPLIFIED  EQUATIONS  FOR  THE  3-D  GAME 


For  most  cases  fuel  consumption  during  the  interception  engagement 
can  be  neglected,  eliminating  Eq.  (3.7).  Assuming  that  thrust  is  nearly 
aligned  with  aircraft  longitudinal  axis  (ftp  «  0)  and  that  angles  of  attack 
are  not  excessive  (implying  cos  a  «  1  and  T  sin  a  «  L)  further 

simplications  are  obtained  in  Eqs.  (3.4)-(3.6). 

i.»o  game  can  be  also  formulated  in  relative  coordinates,  replacing  che  '..iosets 
of  Eqs.  (3  ij-(3  5  j  (for  the  pursuer  and  the  evader)  t>y  only  one  (for  the  relative  co¬ 
ordinates;  However,  only  the  horizontal  components  can  be  meaningful -y  expressed  m 
this  way  since  the  aerodynamic  and  propulsive  forces  actiiig  on  the  airplanes  strongly 
depend  on  their  respective  altitudes.  Consequently,  the  final  set  of  differential 
equations  describing  the  3-D  interception  game  is  of  10  independent  state  variables  . 


(Ax)  =  V£  cos  y£  cos  xE  -  Vp  cos  Yp  cos  Xp 


(Ay)  =  V£  cos  ye  sin  xE  -  Vp  cos  Yp  sin  Xp 


hp  =  Vp  sin  Yp 


h^  =  V,_  sin  Yr 


(3.25) 

(3.26) 

(3.27) 

(3.28) 


np  cos  y 


n  cos  p 


b 


;  =x"Lslnl,P 


P  vp  cos  Yp 


(3 


n„  sin  y_ 

;  ,ii  E 


E  V£  cos  y£ 


(3 


Sometimes  it  is  convenient  to  express  relative  geometry  in  spherical 
coordinates : 


r  =  -Vp  [cos  6  cos  Yp  cos  (Xp-'lO  +  sin  0  sin  Ypj  + 


1 


+  V£  |  cos  0  cos  Yg  cos  +  sin  0  sin  y£ 


6  =  f  {-Vp  [cos  0  sin  Yp  -  cos  Yp  sin  0  cos  [Xp-'J'j 


+  V  cos  0  sin  Yc  -  cos  yc  sin  0  cos 

fc  c  h 


!xe-*) ]} 


<1*  = 


r  cos  0 


-Vp  cos  Yp  sin(xp-4»]  +  VE  cos  y£  sin(x£-V 


(3 


(3 


(3 


In  this  case  an  additional  differential  equation  is  required  for  the 
altitude  of  one  airplane  (hp  for  example),  and  the  opponent's  altitude 
(hE  in  this  case)  is  obtained  from: 

h£  =  hp  +  r  sin  0  (3 


The  relative  geometry  between  the  aircraft  is  depicted  in  Fig  3.1. 


I 


1 

i 

! 


f 


1 


5 .  SIMPLIFIED  MODELS 

From  the  general  model  presented  in  3.4  several  simplified  2-D  models 
can  be  derived. 


5.1  Interception  in  the  Horizontal  Plane 

Motion  in  the  horizontal  plane  requires: 


yp  =  ye  =  e  =  o 


Therefore  the  aircraft  bank  angle  is 


COS  )J  =  — 

n 


Equations  of  the  line  of  sight  variables  in  the  horizontal  plane  are: 


(3.39) 


r  =  -V  cos(xp-<J>)  +  V  cos  (xp-^) 


r 


-VP  sinJxp-’P}  +  VE  sin(xE-^)j 


(3.40) 

(3.41) 


fongitu-Ii  .1  acceleration  of  each  aircraft  is: 


J  *|f 


£T  -  D_  -  n  D , 
max  0  I 


(3.42) 


and  its  turning  rate: 


v  =  £  rn\i'h 

A.  \r  t11 


(3.43) 


The  turning  rate  constraints  are: 

(a)  Below  the  "corner"  velocity  (V  i  V  ) : 


X  <  g 


f  loV  S  C  (M)' 
I  max 


.L 


2W 


2 

-4 

v2! 


(3.44) 
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(b)  Above  the  corner  velocity  (V  *  V  ) 


*  g  r  2 


V  v  max 


This  model  was  used  as  ?  basis  of  the  2-D,  variable  speed  interception 
engagement  reported  in  [3] . 


5.2  Interception  in  the  Vertical  Plane 

In  the  vertical  plane  we  have  the  following  equations  for  the  line  of 
sight  variables. 


r  =  -V  cos(y  -9)  +  V  cos(yp-91 


(3.46) 


9=7  f-vp  sin(yp-e)+  Vg  sin^-B)  j 


(3  4'7’’ 


Longitudinal  acceleration  of  each  aircraft  is: 


r(  T  -  D.  -  n  Dt  i 

- 1  -  sin  t] 


(3-48) 


and  its  rate  of  climb  is: 


Y  =  f-  (n  -  cos  y) 


(3.49) 


The  rate  of  climb  constraints  are  similar  to  (3.44)  and  (3.45)  , 

One  of  the  most  useful  models  for  performance  optimization  in  the 
vertical  plane  is  the  energy-state  model.  In  this  model  the  aircraft  velocity 
is  replaced  by  specific  energy  (E),  defined  by: 


This  transformation  of  variables  leads  to  the  following  equations  for  each 
aircraft. 


ii;  =  V  cos  y 
h  =  V  sin  y 

S  =  (i%-!)v 

Y  =  y  (n  -  cos  Y) 


(3.51) 

(3.52) 

(3  53) 
(3  54) 


Line  of  sight  equations  in  this  model  are  (3.46)  and  (3.47)  n.iilt: 


XP  "  ^  =  XE  "  ^h  =  ° 


(3  55 j 


6.  FORMULATION  OF  THE  DIFFERENTIAL  GAME  AND  THE  NECESSARY  CONDITIONS 
FOR  OFTIMALITY 


As  mentioned  in  the  introduction  we  deal  with  an  interception  problem  in 
which  roles  are  well  defined.  One  airplane  (P)  is  chasing  the  other  (E)  m  a 
3-D  space.  The  problem  can  be  thus  formulated  as  a  zero-sum,  perfect  informa¬ 
tion  game  Assuming  that  capture  is  guaranteed  the  payoff  is  capture  time 

which  P  attempts  to  minimize  and  E  to  maximize. 

Termination  of  the  game  is  defined  by  reaching  a  separation  distance  l 
representing  the  pursuer's  firing  envelope. 


)t2  -  r2(tf)  -  (Ax)2v(Ay)2  ■  fhp-hE)2  j 

! 

i r (tf ) |  <  0  1 


(3.56) 


If  we  assume  the  existence  of  a  saddle  point  solution  the  necessary 


conditions  of  optimality  can  be  applied. 
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First  define  the  Hamiltonian: 


3C  =  1  +  A  x  +  a  y  +  A,  h_  +  Ar  h_  +  Aw  V_  +  A,,  V_  + 
x  y  hp  P  EE  Vp  P  Vp  E 

•  •  •  • 

+  A  >_  +  A  Y„  +  A  v_  +  A  Xt-  +  constraints  (3.57) 

Yp  P  YE  E  xp  P  XE  E 

The  necessary  conditions  of  optimality  include  the  adjoint  equations  and  the 
respective  transversality  conditions. 

The  adjoint  equations  are: 


?  aw 

AX  =  '  37 


{  3W 
xy  =  -  57 


\ '  ^ 


\ =  ■  ^ 


VV  ■  2  Wtf 

(3.58) 

V  V  ■ 

(3.59) 

Xhp(V  *  2thP'hE)tf 

(3.60) 

V(,f)  ■ 

E  f 

(3.61) 

K  (tf)  =  0 

(3.62) 

K  (tf)  -  0 

(3.63) 

t 

A  (t.)  =  0 

YpV  f 

(3.64) 

\  (t f)  -  o 
yE  r 

(3.65) 

x  (t.)  =  0 

Xp  f 

(3.66) 

A  (t.)  =  0 

XR  f 

(3.67) 
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7.  THE  OPTIMAL  STRATEGIES 

Necessary  conditions  for  optimality  are  expressed  by  Isaacs  ME1 : 


min  max  K(5p,EE,yp,yE,np,nE)  =  o  . 


(3.68) 


7 . 1  Pursuer 


The  optimal  throttle  control  for  A  /  o  is: 

VP 


£p  ■  \  i1  ~  sign  \) 


(3.69) 


For  singular  thrust  arcs  (a^  =  L.  =  o)  we  have  0  i  E,*  <  1 .  However, 

P  > 

it  can  be  shown  that  such  singular  arc  is  not  optimal.  Assuming  that: 

T  sin(a+B_,)n  «  Ln,  y  *  is  determined  by: 
maxp  TP  P  ’  P 


tg  b*  = 


A  /cos  Yp 
Xp 


(3.70) 


np  is  determined  by: 


(  2  2  2  i>5  i/cos  Y 

nv  =  Wp  X  +  Av  cos  Yj  - rr- 


P  P  '  Xp  Yp 


P'  2Av  V 
P 


(3.71) 


Equation  (3.71)  is  valid  only  if  the  constraints  on  the  load  factor 
[(3.15)  and  (3.16)]  are  not  violated,  otherwise 


n*  =  inf  n*  ,  n.  ,  n  ' 
P  '  P  Lp  maxp 


(3.72) 
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7.2  Evader 


Similar  expressions  are  obtained  for  the  evader  (maximizer) 


C*  1  (  y  \  } 

SE  =  J{1  *  sign  Ay  ’ 


(3.73) 


tg  y 


A  /cos  y 
Xc  E 

*  E 


(3.74) 


-.ooi-  !/cos  y- 
"E  -  WE  +  \  C0S  YE'“  2U"1T 


(3.75) 


VE  E 


and  if  the  load  factor  is  violated 


ni.  =  inf  r.i£  ,  n,  ,  n 
E  '  c  '  L£  maxE 


(3.76) 


8.  A  CONCLUDING  REMARK 

The  set  of  necessary  conditions  for  game  optimulny  include  20  differential 
equations  (for  the  state  and  adjoint  variables)  and  6  expressions  for  the 
control  strategies,  leading  to  solve  a  nonlinear  two  point  boundary  value 
problem.  For  this  very  complex  problem  the  merits  of  an  analytical  approxima¬ 
tion,  based  on  application  of  a  forced  singular  perturbation  technique,  can  be 
fully  appreciated.  As  a  first  step  of  such  analysis  some  important  modelling 
considerations  are  discussed  in  the  next  chapter. 
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Geometry  of  the  3-D  interception. 


SECTION  IV 

SPT  MODELLING  CONSIDERATION 

The  most  appealing  feature  of  an  SPT  solution  lies  in  its  potential 
to  provide  feedback  control  strategies.  Whenever  such  solution  cannot  be 
obtained  the  SPT  approximation  loses  substantially  from  its  attractiveness. 
Implementation  of  SPT  for  solving  the  multi -dimensional  nonlinear  medium  range 
interception  problems  in  a  feedback  form  depends  on  several  factors. 

(a)  The  ability  to  transform  the  original  problem  to  a  set  of  multiple  time 
scale  (multiple  boundary  layer)  problems  preferably  of  a  single  state  variable 
in  each  layer.  Such  transformation  has  to  be  based  on  observing  an  actual 
time  scale  separation  of  the  original  problem  and  it  requires  a  considerable 
insight  and  sound  engineering  judgement. 

(b)  Terminal  boundary  layers  have  to  be  avoided.  Previous  studies  [14,  15] 
have  indicated  that  if  the  solution  includes  a  terminal  boundary  layer,  this 
part  cannot  be  implemented  in  a  real  feedback  form,  even  if  the  open  loop 
solution  of  the  problem  can  be  computed.  For  an  on-line  implementation, 
intending  to  reproduce  the  open  loop  solution,  the  initial  boundary  layer 
control  strategy  (which  is  uniformly  valid  for  the  boundary  layer  as  well  as 
for  the  free  stream)  has  to  be  switched  to  a  different  terminal  boundary  layer 
strategy  at  a  precise  instant  as  required  by  the  open  loop  matching  condition. 
The  exact  timing  is  critical  in  order  to  satisfy  the  prescribed  end  condi¬ 
tions.  Being  in  the  "free  stream"  no  state  feedback  information  can  be  used 
to  determine  the  conditions  for  the  boundary  layer  initiation  Moreover,  even 
if  it  is  correctly  started  the  terminal  boundary  layer  trajectory  computed  in 
the  real  (forward)  time  direction  is  unstable  with  respect  to  disturbances  in 


the  state  variables.  For  these  reasons  it  is  adviced  to  avoid  formulations 
which  lead  to  a  terminal  ooundary  layer. 

(c)  Use  of  zero  order  approximations  only.  It  was  shown  [15,  16]  that  first 
or  higher  order  correction  terms  improving  the  zero-order  .approximation  have 
to  be  computed  by  integration  along  the  zero -order  trajectory.  It  is  clearly 
an  off-line  process.  Therefore  feedback  implementation  of  SPT  is  inherently 
limited  to  the  zero-order  solution.  These  guidelines  will  be  used  later 
during  the  analysis  of  the  3-D,  variable  speed  interception  model. 

The  extension  of  2-D  constant  velocity  to  3-D  constant  velocity  was 
possible,  because  in  a  constant  speed  model  turning  performance  is  independent 
of  the  position  coordinates.  For  example,  the  extension  of  2-D  constant  velocity  to 
2-D  horizontal  variable  velocity  was  poss  ible  because  (h=const.j  and  V=f(V)  . 

In  a  realistic  3-D  modelling  altitude  plays  a  much  more  important  role  than 
merely  an  additional  position  coordinate  because  of  the  following  reasons: 

I.  Maximum  thrust  is  function  of  both  Mach  number  and  altitude . 

II.  Induced  drag  (caused  by  aircraft  maneuver)  is  inversely  proportional  to 
air  density,  i.e.,  the  capability  to  increase  specific  energy  diminishes 
with  altitude. 

111.  Stall  speed  is  an  altitude  dependent  variable  (which  increases  with 
altitude) . 

IV.  At  low  altitudes  any  attempt  of  an  evading  aircraft  to  escape  a 

descending  pursuer  (attacking  from  higher  altitude) ,  according  to  the 
line  of  sight  strategy  can  lead  to  violation  of  the  ground  state 
constraint  (3.21) . 

V.  At  high  altitudes  any  attempt  to  execute  a  maneuver  with  a  relatively 
high  load  factor  can  force  the  aircraft  (pursuer  or  evader)  to  cross 

the  boundary  of  its  dynamic  flight  envelope. 
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Observing  these  effects  it  seems  that  a  direct  extension  of  the  results  of  a 
- D  i'c:pe3u  modal  to  a  3-D  ono  Lor  al'cerium. v o  1/ adding  thy va.iaolc  speed 

feature  Co  the  3-D  constant  speed  model  as  done  for  the  horizontal  case)  may 
iOt  be  appropriate. 

In  order  to  understand  properly  the  influence  of  the  above  mentioned 
effects  it  is  proposed,  first,  to  analyse  an  interception  engagement  in  the 
vert. cal  plane 

iti  the  past,  energy-state  modelling  has  been  offered  as  an  easily 
Imp! ementable  approach  towards  aircraft  performance  optimization  in  the 
v  : .  ,:;ti  plane.  The  energy-state  approximation  is  based  on  the  observation 
that  in  many  aircraft  maneuvers  the  specific  energy  E  ,  defined  in  (3. 50  , 
varies  slowly  compared  to  the  other  state  variables  [16] .  Thus  the  use  of  a 
reduced  order  model  where  the  "fast"  variables  are  assumed  to  be  pseudo¬ 
controls  seems  to  be  justified. 

This  approximation  has  provided  an  excellent  insight  leading  to  an 
unproved  understanding  of  high  performance  aircraft  trajectory  optimization 
problems.  However,  energy-state  solutions  hav/e  had  only  a  limited  value  for 
direct  an  borne  applications.  As  an  inherent  property  of  the  reduced  order 
model,  discontinuities  of  the  "fast"  state  variables  (like  altitude)  are 
required.  Moreover,  the  model  can  not  satisfy  part  of  the  initial  and 
terminal  conditions  of  the  real  problem. 

A  way  to  overcome  these  inconveniences  is  by  application  of  singular 
perturbation  technique  or  (in  the  absence  of  a  small  parameter  of  physical 
significance)  by  use  of  the  forced  singular  perturbation  technique  (FSPT) . 

In  the  next  chapter: 

fa)  The  suitability  of  the  vertical  plane  medium  range  interception  model 

to  a  FSPT  formulation  will  be  checked. 

(b)  The  model  will  be  solved  in  detail. 
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SECTION  V 


FSPT  SOLUTION  FOR  MEDIUM  RANGE  INTERCEPTION  IN  THE  VERTICAL  PLANE 


1.  GENERAL  CONSIDERATIONS  AND  EQUATIONS  OF  MOTION 

Investigating  the  suitability  of  this  problem  to  a  FSPT  formulation,  it 
is  first  necessary  to  identify  "fast"  and  "slow"  variables  in  the  equations 
of  motion-  For  this  purpose  we  assume  the  following: 

(a)  Medium  range  interception  is  characterized  by  a  separation  distance 

largely  exceeding  the  aircraft  turning  radii.  As  a  result,  the  rate  of 
change  of  the  relative  geometry  (expressed  by  the  line  of  sight  components) 
is  the  slowest  variable  in  the  model. 

1  (b)  dased  on  engineering  observations,  specific  energy  is  slow  variable 

compared  with  V  or  h  [17].  Therefore,  use  of  E,h  instead  of 
V,h  is  essential  towards  legitimate  use  of  FSPT. 

,  (c)  The  fastest  variable  in  this  model  is  the  climb  angle  y  (since  h  is  an 

integral  of  y) . 

Applying  these  considerations  to  the  vertical  plane  interception  (whose 
geometry  is  depicted  in  Fig.  5.1),  the  equations  of  motion  (with  £  inserted 
j  artificially)  are: 


4 


4 


4 
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(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 


1*1,  <  0  (5.9) 

f 

where:  l  -  Pursuer's  capture  range. 

The  natural  payoff  is  time  of  capture  which  the  pursuer  tries  to  minimize 
and  the  evadeT  to  maximize. 

Before  starting  the  formal  solution  several  additional  assumptions  ought  to 
be  made. 
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(a)  Interception  time  is  long  enough  to  allow  each  aircraft  to  reach  its 
optimum  altitude  (at  which  maximum  velocity  can  be  achieved).  Later  on 

it  will  be  shown  that  this  phase  corresponds  to  the  reduced  game  solution. 

(b)  Pursuer's  capture  range  is  greater  than  its  turning  radius.  Therefore 
the  "end  game"  phase  disappears  and  the  interception  terminates  with  a 
"tail  chase"  configuration. 


.2?  _  W 

hP  Ep=const . 


1  Ahp(tf)‘  2(VhE)tf  (5‘14^ 


W 

E  E  =const. 
b 


+  *  AhE(tf)=-2(hp-hE)tf  (5-15) 


;  \  (tf)  =  o 

rp  r 


(5.16) 


•E  ' 


;  A  (tf)  -  o 


(5.17) 


Optimal  strategies  are  derived  from  Isaacs  ME1 


min  max  ?f  =  0 


^p*np  ^e ’^E 


2.1  Pursuer 


Sp  =  j  (1  "  sign  Ae  -1 


A_  ^  0 
P 


(5.18) 


It  can  be  shown  that  singular  arcs  (A  =  A  =0)  are  not  optimal. 


*  _  oo\ 

is  calculated  from  —  =  0. 


Ip  Wpg 


p  \  ■>*  v2 

Lp  “l IPVP 


(5.19) 


Equation  (5.19)  is  valid  only  if  the  constraints  of  the  load  factor  are  not 
violated.  A  more  general  expression  for  n*  is: 


nP  '  lnf(nP-%'\axp) 


(5.20) 
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2.2  Evader 


*  J  C1  ♦  sign  XE  ) 


Xp  *  0 
fcE 


(5.21) 


Similar  to  the  pursuer,  singular  thrust  arcs  are  not  optimal  in  this  case. 

nl  is  calculated  from  —  =  0. 

E  on_ 

E 


*  ^E  V 

"E  =  XE  2D  V2 
E  I  E 
E 


(5.21a) 


or,  when  constraints  are  taken  into  consideration: 


nE  ’  n  ) 

E  E 


(5.22) 


3.  THE  REDUCED  GAME 

By  setting  e  =  0  in  equations  (5.2)  f  (5.7)  and  (5.11)  t  (5.17)  the  reduced 
game  is  obtained.  Denoting  all  variables  of  the  reduced  game  by  the  superscript 
"o"  we  get  from  (5.2)  '■  (5.7): 


-°  T 
"P  maxr 


f  °  T 
f:  maxr 


°2  o 

"p  Di. 


o2  o 
nE  PT 


(5.23) 


,,0  .  0 
vp  Sln  Yp 


1,0  •  .  0 
Vp  sin  Ye 


y°E  -  o 


(5.24) 


4 


g  ,  o  o- 

£  Cnp  -  cos  yr, 


"p  -  COS  Yp 


(5.25) 


-i  (n°  -  cos  v°) 
VE 


o  o 

nc  =  cos  yc 
E  E 


"5e~  =  "jg- 

'P|hp=const.  E  h^const. 


(5.26) 


-^7—  =  'TT — 

P  Ep=const.  nE  E^=const. 


(5.27) 


3jK°  djf 


(5.28) 


The  controls  in  this  game  are  the  specific  energy  and  altitude  of  each 
aircraft.  These  controls  are  obtained  by  min-max  of  the  reduced  game  Hamiltonian. 


*°  =  1  +  (-v°  ♦  V°)  =  0 

h 


(5.29) 


From  (5.11)  and  (5.9)  it  follows  that  >°  >  0.  Using  this  fact  and 

rh 

performing  maximization-minimization  of  the  Hamiltonian: 


hp,F°  =  arg  min  {iH0} 

VEP 


(5.30) 


hp,E°  =  arg  max  {/} 


(5.31) 


hE,EE 


we  obtain: 


=  arg  max|v°} 


(5.32) 


hr»Ep 
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(5.33) 


h“,E°  =  arg  max{v°} 

Iv.E,, 


In  other  words,  to  find  the  altitude  and  energy  level  at  which  each 
aircraft  flies  during  the  reduced  game  one  has  to  search  the  whol;?  range  of 
energies  and  altitudes  of  each  aircraft's  flight  envelope. 

Since  Yp  =  Yp  =  0  the  reduced  game  is  characterized  by  flight  at  constant 
altitude  at  which  maximum  velocity  can  be  achieved.  This  velocity  is  determined 
by: 


V°  =  sup  arg  {£p  T 


P  max 


sup  arg  T 


max 


From  which  we  obtain:  = 


1 


V°  ♦  V° 
P  E 


-  D°  =  0} 

P 

-  D°  -  0} 

E 

And  finally 


(5.34) 

(5.35) 


(5.36) 


4.  THE  OUTER  BOUNDARY  LAYER  (ENERGY  BOUNDARY  LAYER) 

The  equations  of  motion  in  the  outer  boundary  layers  (separately  controlled 
by  each  player)  are  obtained  through  the  stretching  transformation:  =  t/£. 

Denoting  all  variables  in  the  outer  boundary  layer  by  the  superscript  "i" 
we  have,  setting  £  =  0  and  using  the  matching  conditions: 


(5.37) 


maxr 


"p2  v 


(T^O) 


(5.38) 
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dEE  VE  i  i 

_ E  =  _E  (F  T 

dTi  *;  E  mxe 

E 


D»E  -  V 


i  .  i  i 
0  =  £  dFT  =  VP  sin  yp  ~  YP 

i 


e  *r  ’  ve  si"  ye  -  ye 

x 


0  - £2  ar  -  4  ("p  - cos  ~  *1 

1  vp 


dYc 


0  =  e 


2  "* 1  E  K  ,  l  _.i*  i 

=  -4"  (n_  -  cos  Yc)  ~  nc 


dx. 

l 


=  1 


V. 


d  r 

il  =  0  -  A1  (r.) 
dx.  r  l 

l  h 


=  A 


3 

dx. 

l 


_  _g_ 

vp  3vp 


K 


dt. 

1 


11  =  X 

v1  av1 

v  E  E 


0  =  e 


0  =  E 


J. 

hp  _ 

ay1 

dx . 

i 

ahp! 

dA  i 

E  .. 

a^ 

dx. 

i 

ahj 

E 

Ep=const 


E  =const , 
E 


(T.=0) 


(5.39) 

(5.40) 

(5.41) 

(5.42) 

(5.43) 

(5.44) 

(5.45) 

(5.46) 

(5.47) 

(5.48) 


31 


(5.49) 


0  =  e 


0  =  e 


dX  i 
2  P 

dx. 

l 


dX  i 

2  _ 

dx. 

l 


93^ 

*p 

331* 

*E 


(5.50) 


In  this  boundary  layer  we  still  have: 


Each  aircraft  has  to  change  its  energy  level  from  its  initial  value  (E  or 

0 

Ec  )  to  the  value  dictated  by  the  outer  solution  (E°  or  E°) .  To  do  it 
fc0  F  b 

optimally,  it  is  necessary  that  at  each  energy  level  the  aircraft  will  be  at 

an  optimal  altitude  (considered  as  control  in  this  boundary  layer) .  This 

altitude  is  calculated  (for  both  aircraft),  by  minimization-maximization  of  the 

Hamiltonian  GH1)  . 


V1  =  1  +  X°  (-v£  ♦  vi)  +  X* 

\  bP  bP 


*  \l  Pp 
nj=l  EE  SE 


i  , 
n  =1 
E 


=  0, 


(5.51) 


A  • 

where  P  =  E  is  the  specific  power. 

Substituting  X°  (5.36)  into  (5.51)  we  obtain: 


V°  -  V1 

ffl  P  P  l1  n1 

=  — -  — ~  +  X  P 

v?  -  v°  Ep  sp 

r  £ 


v°  -  V1  .  , 

•  -J — 

"i-i  v?  -  ve  e  e 


nE=1 


=  0 


(5.52) 


Since  the  Hamiltonian  is  separable  (with  respect  to  the  players  controls) : 


V1  -  V° 
P  P 


1 


ep  v°  -  v°  pi 

P  E  Sr 


<  0 


(5.53) 
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-v1  +  v° 

E  E  1 

V°  -  V°  p1 

P  E  S 


>  0 


The  optimal  altitudes  (hp,h*)  are: 


hp  =  arg  min  W1 
'  hP 


hP  =  arg  max 
hE 


}& 


Expressed  more  specifically. 


h  =  arg  max 
hP 


h  =  arg  max 
“  hE 


From  (5.53),  (5.54): 


=  £*  =  1  . 


(5.54) 


(5.55) 


(5.56) 


(5.57) 


(5.58) 


(5.59) 


5.  THE  INTERMEDIATE  BOUNDARY  LAYER  (ALTITUDE  BOUNDARY  LAYER) 

The  equations  of  motion  in  the  intermediate  boundary  layer  (separately 

controlled  by  the  players)  are  obtained  through  the  stretching  transformation 
2 

i j  =  t/£  .  Denoting  all  variables  in  this  boundary  layer  by  the  superscript 
"j”  we  have,  setting  e  =  0  and  using  the  matching  conditions: 


dr J 

h  l 

dT .  u  rh 
J 


(5.60) 
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EjUyoj 


EE(Tr0) 


dh,, 

d?”  “  Vl  sin  YP 


VE  Sl"  YE 


-W 

\/J  1 


dvJ 

(nJ 

drj  V  E 


hjtyo,  =  h 


hE^r0)  =  \ 


cos  Yp) 


cos  ->£) 


nJp  =  cos  yJp 


nJE  -  cos  Yg 


dXr 

rh 

dT. 

J 

=  0 

rh 

(V 

=  i 

dxj 

fcP 

dT 

J 

=  0 

‘V 

—  / 

dA^ 

lh 

dr. 

J 

=  0 

"  % 

(V 

=  / 

d\l 

hP 

AJf1  + 

X 

930* 

d'l  . 

J 

3hj 

vj 

p 

9vj 

(5.61) 


(5.62) 


(5.63) 


(5.64) 


(5.65) 


(5.66) 


(5.67) 


(5.68) 


(5.69) 


(5.70) 


i 


(5.71) 


dtf  t  g  dJfJ 

3vi 


0  =  e 


(5.72) 


0  =  F. 


(5.73) 


Velocities  in  this  boundary  layer  are  calculated  by: 


Yj  {2g(Ej  -  hj)} 


(5.74) 


According  to  (5,72),  (5.73)  the  path  angles  (Yp,Yg)  become  the  controls  of 
this  boundary  layer.  Their  role  is  to  transfer  the  airplanes  from  their  current 
state  (hp,hg)  to  the  required  "optimal"  altitudes  (hp,h*)  determined  by 
(5.57),  (5.58). 

The  Hamiltonian  of  this  boundary  layer  is: 

^  =  1  *  x°  (-vi  cos  yj  ♦  vj.  cos  ‘ril  *  ij-  (Tji,-pb  •  X*  ^  (Tg-Dg)  * 

%  V  V3 


*  Xl  Vp  sin  YP  *  xh_  VE  Si"  YE  ■  °' 

i  E 


(5.75) 


Using  (5.72),  (5.73)  we  have: 


A°  vj  sin  yj  +  AjJ  vj  cos  yJ  =  0 
ho  p 


(5.76) 


-a  vl  sin  'A  +  cos  y^ 

r,  fc  E  hn  E  t 

h E 


(5.77) 


After  using  (5.29),  we  obtain: 
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(5.78) 


xj  =  ~l8--p- 

hF  v°-v° 


tg  ye 


hE  V°  -  V° 

C  VP  E 


A->  = 


(5.79) 


Substituting  A°  (5.29),  A*  (5.53),  A*  (5.54),  A^  (5.78),  Ajj  (5.79)  into 
Th0  %  %  P  E 

the  Hamiltonian,  noting  that  3C-1  =  o  and  using  the  fact  that  the  Hamiltonian  is 

separable  into  two  parts  (exclusively  controlled  by  the  pursuer  or  the  evader) 


we  have: 


Vj  P1  1 

P  SP 

Y„  =  cos  1  - = - t* — - * - s-  'sign  (h^  -  hjh 

P  V^P1  -  )  ♦  V1  P^  P  p 

PL  sp  spJ  VPQ  Ksp 


(5.80) 


yj  p* 

E  KSr 


Yp  =  cos 


">1  - *  vL  i 


sign  (hg  -  hj) 


(5.81) 


It  has  been  shown  in  | 1 8 J  that  the  arguments  in  equations  (5.80)  and  (5.81)  are 
bounded  between  zero  and  one. 


.6.  THE  INNER  BOUNDARY  LAYER  (TURNING  BOUNDARY  LAYER) 

The  equations  of  motion  in  the  inner  boundary  layer  (separately  controlled  by 

3 

each  player)  are  obtained  through  the  stretching  transformation  =  t/e  . 
Denoting  all  variables  in  this  boundary  layer  by  the  superscript  "k"  we  have, 
setting  e  =  0  and  using  the  matching  conditions: 


ur,  , 


(5.82) 
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dK 

-3—Z-  =  0  -  A*  (T,  )  =  A, 

dTk  V  k  ? 


(5.93) 


1 


h 


b 


d\ 


dT, 


,,k 

dA 


a  i , 
K 


3?^ 

3^ 

3Y, 


(5.94) 


(5.95) 


V  ■» 

Velocities  in  this  boundary  layer  are  given  by:  V,  ^  =  VJ  =  V 

w  u 

The  Hamiltonian  of  this  boundary  layer  is: 


3<k  ’  1  4  X“h(-VP0  C°S  VP  *  \  C°S  Yh)  *  %g  V  (Vax„-  D0d-  "p  D  Jo  * 

VE  /  2 

AEc  \Tmaxp"  V  "E  %)<>  +  ^  V 

Yp)  +  X^r  vf~  (nE  -  C0S  Ye)  ’  t5'96) 


~P 


k  xj  V 


'e  E 

fco 


P  Sin  YP  +  Ah  VE  Sin  YE  + 

P0  po 


,  k  g  (  k 

+  \  v  Vnp  '  cos 
p  p  x 
0 


'E  E, 


k  k 

In  this  boundary  layer,  Op  and  n^  are  the  controls.  Expressions  for  their 
optimal  values  are  obtained  as  follows: 


((n. 


TT 

9nc 

E 


(5.97) 


From  which : 


gWr 


nP  = 


2DI  VP 
P0  ° 


(5.98) 
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r 

L 

L 


♦ 


I 


xk 

k  ye 


gwc 


= 


E  Aj  2Dt  Vp 

ee  *e  eo 
L0  0  u 


(5.99) 


Substituting:  A  (5.36),  A  (5.53),  A*  (5.54),  A^  (5.78),  A^  (5.79) 

rh  Cp  Kc  "d  bp 

o  0  po  bn 

into  the  Hamiltonian,  using  its  separability  (with  respect  to  each  aircraft 
k  k 

controls)  A  and  A  can  be  extracted  from  it  and  substituted  into 

YP  YE 

(5,98)  and  (5.99). 


n?  =  cos  Yp  +  i 


0 


cos  yJ  (v°  -vj  )  DI 
0  0  0  1 


[l  -  cos  (yj  -Yp)] 


2  k  2  i 

+  cos  Yp  -  cos  Yp 


\  sign  (Yp  - Yp) 
0 


(5.100 


r 


v 


k  k 

n*  =  cos  «*  - , - - - : —  — 5- — 

1  E  cos  4  K -4  ) 

Lo  fco  fco  h 


■  fl  -  cos  -Yk)] 


0 


2  k  2  j 

-1  COS  Y_  -  COS  Yr- 
r.  E 


0 


sign  (Yp  -  Yp) 

Eo  E 


(5.101 


7.  FEEDBACK  SOLUTION  OF  THE  PROBLEM 

iji.ifoniily  valid  feedback  expressions  for  np  and  n£  can  be  obtained 
if  we  replace  the  initial  (stationary)  values  of  the  various  variables  in 
(5.100),  (5.101)  by  their  current  (measurable)  values. 

The  pursuer's  aerodynamic  load  factor  will  be  expressed  as  follows: 
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np  =  c°s  Yp 


!  V*  Tmaxp'  DOp  '  DIp^ 1  r 

cos  Yp  V°  -  Vp  DIp 


+  COS2  Yp  -  COS2  Yp  ■  *  sign  (Yp  -Yp) 


cos  (Yp  -  Yp)  I  + 


(5.102) 


,9  ?  i.. 

The  term  (cos  Yp  -  cos  YpJ  is  small  compared  with  the  first  one  under  the 
root  sign  as  long  as  Yp  ^  Yp*  When  Yp  -*•  Yp  the  solution  has  to  approach 
asymptotically  a  straight  line  path  at  Yp*  These  two  arguments  justify 


:ancellation  of  this  term. 


We  get  finally: 


,  1  p 

n  =  cos  Yp  +  - J  “o - T 

C0SYp  Vp  -  Vp 


*  («r  W  V  '  "i 


P  P  P 
D, 


1  -  cos  (Yp  -  Yp) 


sign  (Yp  -  Yp) 


(5.103) 


where : 


Vp  -  pursuer's  maximum  velocity  in  a  straight  level  flight  within  its 
flight  envelope. 


Vp  =  arg  {T  -  D  -  D  =  0} 

tp,hp  maxp  °P  lP 


(5.104) 


The  superscript  "i"  is  used  to  denote  variables  obtained  at  h1  (and  the 
current  specific  energy  level)  defined  by: 


Pi  (E) 

i  SP 

h  (H)  -  arg  max  — - i - 

hp  V°  -  Vp(E) 


(5.105) 


P~  (H)  -  Pursuer's  specific  power  at  straight  level  flight  at  an  altitude  hp. 
bP 
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vi  is  calculated  from: 


Vj  =  { 2g  CEp  -  hp) } 


(5.106) 


The  flight  path  angle  Yp  at  which  the  aircraft  has  to  fly  in  order  to 
reach  the  altitude  h1  is  expressed  in  feedback  form  as  follows: 

signfh1  -  h)  (5.107) 


i  -i 
fp  =  cos 


v1  p1 

P  Sr 


v°p(psp-psp)  *  VP  psp 


P  -  Pursuer's  specific  power. 

SP 

iTie  pursuer's  optimal  throttle  remain  as  in  (5.59): 


'P  =  1  • 


Similar  results  hold  for  the  evader. 


8.  TERMINAL  PHASE  OF  THE  VERTICAL  INTERCEPTION  MODEL 

According  to  the  solution  presented  in  this  chapter  the  interception 

terminates  at  the  free  stream  conditions,  i.e.,  both  airplanes  fly  horizontally 

(Y  =  0)  at  their  maximum  speed  at  their  respective  optimal  altitudes  (h°,  h°)  . 

If  h°  t  h°  (which  is  generally  true)  the  transversality  conditions  on 
P  r. 

X,  (t  )  and  X  (t  )  (5.14),  (5.15)  are  not  satisfied  by  (5.78)  and  (5.79). 

hp  t  FE  t 

Relating  to  practical  terms,  if  the  altitude  gap  | hp  -  hE(  is  too  large, 

weapon  firing  may  not  be  convenient.  Moreover,  since  the  choice  of  optimal 

altitude  h°  is  independent  of  the  position  of  the  other  airplane,  such 

strategy  cannot  be  considered  feedback  optimal  for  the  pursuer.  Even  if 

] Ah° (  •'<  t,  the  evader  could  avoid  capture  by  flying  at  a  non-optimal  altitude 

out  of  the  firing  envelop.  For  all  these  reasons  the  energy  state  FSPT  model 
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requires  a  terminal  boundary  layer.  As  it  was  mentioned  previously  (Ch.  4) 
feedback  control  of  a  terminal  boundary  layer  is  not  feasible.  This  difficulty 
can  be  overcome  by  using  a  dual  modelling  approach  for  the  pursuer  only.  The 
evader  has  to  fly  always  at  its  optimal  altitude. 

If  altitude  influence  could  have  been  neglected  (which  is  unfortunately 
not  true),  the  pursuer's  strategy  would  have  been,  similarly  to  the  horizontal 
case,  a  line  of  sight  pursuit.  It  seems  therefore  reasonable  to  propose  for 
the  pursuer  a  weighted  combination  of  these  control  strategies  obtained  from 
'.wo  different  FSPT  models.  This  idea  is  presented  mathematically  by  the 
following  expression  for  the  pursuer's  aerodynamic  load  factor. 


nP  =  Bn  +  (l-B)n 
V 


(5.108) 


where: 


B  = 


(5.109) 


rip  -  Pursuer's  vertical  plane  load  factor  derived  from  line  of  sight  model; 

;ip  -  Pursuer's  vertical  plane  load  factor  derived  from  energy-state  model; 
n  -  Empirical  constant. 

The  time  varying  weighting  factor  8  ensures  a  smooth  transition  and  enforces 
the  condition  of  capture  (i.e.,  convenient  weapon  firing).  It  is  easy  to  see 
that  when  r  =  r^  (at  the  beginning  of  the  interception) 

8  =  0  -  n  =  n 
V  f 

i.e.,  energy-state  model  is  used  exclusively.  When  r  =  r^.  (at  termination) 

8  =  1  -  n  =  n  , 

rv  r 
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i.e.,  line-of-sight  model  is  used  exclusively.  At  intermediate  values  of  r 
(r  <  r  <  r  )  n  is  determined  by  (5.108). 

U  Tty 

fhe  value  of  m  -  2  was  found  acceptable  based  on  some  numerical  experimen¬ 
tation.  In  Figure  5.2  the  final  time  vs .  m  is  depicted  (for  the  example  that 
will  be  presented  in  detail  in  Chapter  7)  showing  that  for  m  -  2  the  pursuer 
minimizes  the  capture  time. 

It  can  be  summarized  that  the  solution  of  the  vertical  interception  game  is 
a  pair  of  uniformly  valid  strategies  expressed  differently  for  the  pursuer 
and  For  the  evader.  These  results  will  be  interpreted  in  the  next  chapter  in 
relation  to  the  3-D  variable  speed  medium  range  interception  game. 


43 


Fig.  5.1:  Geometry  of  the  vertical  plane  interception. 


Fig.  5.2.  Final  time  as  function  of  n. 


SECTION  VI 


[ 

S 


t 


i 


FSPT  STRATEGIES  FOR  THE  3-D  INTERCEPTION 

1.  REVIEW  OF  THE  MODELLING  CONSIDERATIONS 

As  explained  in  Chapter  4,  the  implementation  of  the  methoo  of  forced 
singular  perturbations  to  solve  the  original  3-D  variable  speed  problem 
can  by  no  means  be  considered  as  a  straightforward  extension  of  previous 
results.  The  relative  complexity  of  the  FSPT  game  analysis  for  the  vertical 
interception  presented  in  Chapter  5,  compared  to  the  horizontal  case  [2] 
clearly  indicates  that  a  much  more  careful  modelling  approach  may  be  required 
For  an  enhaced  insight  it  is  useful  to  recall  the  models  used  in  the  2-D 
(horizontal  and  vertical)  problems. 


TABLE  1:  FSPT  Models  for  Horizontal  and  Vertical  Interception 


Horizontal  Interception 
(Model  A) 


Vertical  Interception 


|  r  -  rh  -  -Vp  cos  (xp V£  cos(xE-^) 
$  =  “  [-Vp  sin(xp-U»)  ♦  VE  sin(xE-^) 

*  /Tp_r>pN 

> 

•  /tf"de\ 

cv*  -hr)* 


P2v  ,  g  (n2 

e  xp  *  vj^np  '  i j 


s  X  fn2  .  i]*5 
V  v  E  J 
E 


rh  =  -vi>  cos  Yr  *  ve  cos  ye 


fEE  *  (tt)  VE 


c~hp  =  Vp  sin  Yp 


£thE  =  VE  Sin  XE 


3*  g  r 
e  Y~  =  rf*-  n 


p  -  vj  -"P  •  C0s  V 


e^YE  -  £  ("e  ■  C0S  YE! 

M 
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Inspection  of  Table  1  suggests  the  following  hierarchy  for  the  3-D  model. 

a)  Horizontal  components  of  the  relative  geometry  (Ax, Ay  or  r^,<Jj)  are  no 
doubt  the  variables  of  the  reduced  game  for  a  medium  range  interception. 

b)  Since  in  the  horizontal  model  (model  A)  speed  can  be  expressed  by  specific 
energy,  there  is  no  contradiction  between  the  two  models  to  appoint  the 
specific  energy  as  the  first  (slowest)  boundary  layer  variable. 

c)  The  experience  with  energy  state  models  certainly  requires  to  admit  that 
altitudes  vary  faster  than  specific  energy. 

d)  The  flight  path  angle  is  a  definitely  faster  variable  than  altitude  (since 
h  is  the  integral  of  y) . 

e)  Turning  dynamics  in  both  planes  seem  to  belong  to  the  same  level  of 

hierarchy,  being  the  components  of  the  same  variable:  the  aircraft  total 

— ► 

angular  velocity  vector  w  given  by 
2  •->  2 

go  =  y“  +  (x  cos  y)  (6.1) 


Substituting  to  (6.1): 

y  =  f  (ny  -  cos  y) 


A 

ny  =  n  cos  u 


(6.2) 


and 

‘  ^ 

x  V  cos  y 


A 

n^  =  n  sin  y 


(6.3) 


leads  to: 


2n  cos  y  cos  y  +  cos  y 


(6.4) 


Taking  into  account  that  the  total  aerodynamic  load  factor  "n"  is  a 
constrained  control  variable  by  (3.15)  and  (3.16)  the  coupling  between  the 
horizontal  and  vertical  components  must  not  be  overlooked. 
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However,  for  unconstrained  load  factor  the  exact  optimal  solution  (i.e.. 


Eqs . (3. 70) - (3. 71)  can  be  expressed  by  defining 


(6.5) 

(6.6) 

(6.7) 

(6.8) 


These  expressions  suggest  that,  for  cos  Y  w  1,  the  vectorial  composition 
of  the  two  perpendicular  2-D  solution^  can  be  a  reasonable  approximation  for  a 
3-D  problem  subject  to  properly  imposed  cons  raints. 

Based  on  this  insight  it  is  therefore  proposed  that  the  zero-order  FSPT 
solution  of  the  three-dimensional  variable  speed  medium  range  interception 
game  will  be  synthesized  by  this  approach. 


2,  OPTIMAL  CONTROL  STRATEGIES  FOR "3-D  MEDIUM  RANGE  INTERCEPTION  _ 

The  zero-order  composite  FSPT  strategies  of  the  aircraft  can  be  now 
summarized. 

For  both  players  the  optimal  throttle  setting  is  always  maximal. 


( rC''  _  CrC'l 

'  V  3D  "  3D 


1 


(6.9) 


Let  us  note  in  this  conjunct  ion  that  a  model  B  [3]  type  solution  is  not 
necessary  in  the  3-D  case.  Airplanes  can  loose  speed  rapidly  (if  it  is 
required)  by  exchanging  kinetic  energy  to  altitude  without  giving  up  total 
specific  energy. 
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The  aerodynamic  load  factor  and  the  respective  bank  angle  are  given  by 


,  c*  .  f  ,2  2  \ 

V3D  ‘  "ln  \  *  J 


(6  10) 


tg  yp  =  /tv 
r  p  p 


(6.11) 


with  similar  expressions  for  the  evader.  The  horizontal  load  factor  expression 
is  the  same  for  both  players  as  derived  in  [3] 


„  T  -  D.  -  D  1*5 

V  max  0  x  r ,  (  .  1 1  ■  r  i  i 

_ —  .  - - -  [1  .  cos  jM<,  If  sign'^-x. 

V  -V  DI  J 


(6.12) 


It  has  been  however  remembered  that  the  vertical  maneuvering  strategy  of 

the  pursuer  is  different  from  the  evader's,  riy  is  given  by  Eqs .  (5  109)- 

VP 

(5.110),  while  ny  is  simply  n£  expressed  by  Eq.(5.103). 


SECTION  VII 


A  DEMONSTRATIVE  EXAMPLE 


1.  AIR  DEFENSE  SCENARIO 

In  order  to  demonstrate  the  applicability  of  the  zero-order  feedback  solution 
as  outlined  in  Chapter  6  an  imaginary  air  defense  scenario  was  selected. 

Aircraft  "A"  (the  pursuer  in  the  game),  anagile  air  superiority  fighter  of 
high  thrust  to  weight  ratio,  patrols  to  defend  a  small  area  from  an  eventual  attack 
of  a  tactical  bomber  "B"  (the  evader)  penetrating  at  a  low  altitude  at  high  sub¬ 
sonic  speed. 

Some  basic  information  on  both  aircraft  is  given  in  Table  7.1,  Aerodynamic 
and  thrust  data  of  aircraft  "A"  was  obtained  and  used  in  the  computer  program  in 
a  tabular  form  as  shown  in  Table  7.2.  The  data  of  aircraft  "B"  is  expressed  by 
a  polynomial  approximation  in  the  computer  program  and  depicted  in  Figs.  7.1  and 
7.4. 

Comparison  of  the  flight  envelopes  for  both  aircraft  (Fig.  7.5)  shows  a  definite 
advantage  of  "A"  in  low  altitude.  "B"  has  however  a  higher  maximum  speed  and 
maximum  altitude. 

It  is  assumed  that  the  interception  game  starts  by  simultaneous  detection 
followed  by  an  immediate  bomb  load  release  by  the  penetrating  "B".  Since  "A"  is 
equipped  with  a  superior  weapon  system,  "B"  wishes  to  enhance  the  survival 
probability  by  optimal  evasion.  The  interception  is  successful  if  the  pursuer 
can  reach  the  firing  range  of  its  missiles  in  finite  time.  If  such  defined 
"capture"  can  take  place, the  objective  of  the  pursuer  is  to  minimize  the  time  of 
capture  and  the  evader  wants  the  maximize  it.  Due  to  the  differences  in  air¬ 
craft  flight  envelopes  (see  Fig.  7.5),  capture  is  not  always  guaranteed. 
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This  scenario  serves  a.-  hast-  for  the  following  studies, 
a  Detailed  time  history  of  a  characteristic  single  engagement  analyzing 
aircraft  strategies. 

b.  Parametric  study  of  the  influence  of  initial  conditions,  weapon  and 
airplane  characteristics  on  the  outcome  of  the  interception. 

Such  a  systematical  investigation,  requiring  to  solve  a  large  number 
of  "optimal"  engagements,  could  be  performed  with  the  limited  computational 
resources  only  due  to  the  inherent  efficiency  of  the  feedback  solution. 

2.  TIME  HISTORY  OF  A  CHARACTERISTIC  ENGAGEMENT 

The  initial  and  end  conditions  of  this  engagement  are  summarized  in 
Table  7,3  and  the  projection  of  the  initial  geometry  are  depicted  in  Fig. 

7.6  In  figures  7.7  -  7.13  the  time  histories  of  the  control  and  state 
variables  are  depicted.  Based  on  these  results  several  observations  can  be 
I  made . 

a.i  During  the  initial  phase  of  the  engagement,  the  pursuer  performs  a 
sharp  diving  turn  increasing  its  velocity  and  closing  the  angular  difference 
with  the  line  of  sight.  At  the  same  time  the  evader  engages  in  a  sharp 
climbing  turn  at  almost  constant  velocity. 

b)  About  35  seconds  after  the  maneuver  started,  both  azimuth  angles  xp 

and  approach  the  line  of  sight  azimuth  angle  and  the  aircraft 

maneuvers  in  the  horizontal  plane  are  almost  completed. 

c)  From  that  point  till  the  end  of  the  interception  the  aircraft  fly,  only 
with  the  intention  to  increase  their  velocity. 
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d)  At  about  t  =  60  sec,  the  evader  starts  diving  in  order  to  change  its 
wiocity  from  subsonic  to  supersonic,  while  at  the  same  time  the  pursuer, 
flying  already  at  supersonic  velocity,  climbs  toward  the  evader.  This 
maneuver  "seems"  to  be  non  reasonable,  however,  according  with  energy-state 
model  this  is  the  best  way  for  the  evader,  to  increase  its  specific  energy 
and  velocity,  required  to  escape. 

e)  It  should  be  noticed  that  the  pursuer  has  more  favorable  initial  condi¬ 
tions  than  the  evader,  since  its  initial  dive  maneuver  serves  two  goals  at 
the  same  time,  increases  its  velocity  and  biing~  -ioser  to  the  evader. 

f)  Another  interesting  point  to  be  noticed  is  that  although  the  maneuver 
starts  with  the  evader  being  at  5000  meters  lower  than  the  pursuer,  at 

t  =  32  sec  both  aircraft  are  at  the  same  altitude,  3000  meters,  and  from 
that  moment  until  termination  of  the  interception  the  evader  flies  at  higher 
altitude  than  the  pursuer. 

g)  As  one  can  see,  capture  occurred  in  this  example,  before  any  one  of 

t..  . ircraic  reached  the  "outer"  (steady  state)  conditions,  i.e.,  V°  = V 

max 

Since  the  strategies  of  both  aircraft  are  based  on  the  assumption  of 
reaching  maximum  velocity  during  the  reduced  game,  and  these  maximum  velocities 
are  not  achieved, the  accuracy  of  the  FSPT  approximation  might  be  degraded 
This  will  be  discussed  in  Chapter  10. 

Observing  the  time  histories  the  following  conclusions  can  be  made: 

1  The  results  confirm  the  assumptions  of  the  FSPT  modelling  about  time 
scale  separation 

2.  After  the  horizontal  turn  is  completed  the  behaviour  of  both  airplanes 
are  dominated  by  the  energy  state  model  requirements.  This  fact  can  be 
clearly  seen  in  Figs.  7  14,  7.15  where  the  optimal  "energy  climb"  profiles 
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h 1 C E)  and  the  actual  trajectories  are  plotted  into  the  flight  envelope  maps 
3  Due  to  the  small  altitude  difference  between  the  airplanes  the  line  of 
sight  strategy  used  by  the  pursuer  at  the  terminal  phase  cannot  be  noticed 


3 .  PARAMETRIC  INVESTIGATION 

The  qualitative  and  quantitative  outcome  of  the  game,  i  e., 

(i)  can  capture  be  achieved? 

(ii)  the  value  of  the  "time  of  capture", 
depend  on  the  parameters  of  the  airplanes  and  the  weapon  system  as  well  as  on 
the  initial  conditions.  For  a  brief  but  systematical  analysis  the  influence 
of  the  following  were  considered. 

First  the  influence  of  initial  range,  pursuer’s  initial  heading  and  missile 
firing  range  (capture  radius)  on  the  qualitative  outcome  was  investigated. 

In  Fig  7  16  the  relevant  capture  zone  projections  for  different  firing  ranges 
are  depicted.  It  shows  that,  for  the  set  of  initial  conditions  kept  constant, 
capture  can  be  attained  even  in  the  most  unfavourable  conditions  for  the 
pursuer  (xp  =  18()u  )  . 

Figs.  7.17,  7. Id  show  the  influence  of  the  initial  altitude  on  the  time 
of  capture,  indicating  that  almost  independently  of  initial  range  and  capture 
radius  the  best  altitude  for  patrolling  is  between  7-8  km  This  gives  to  the 
pursuer  a  2:1  advantage  in  specific  energy.  The  results  show  that  if 

h  <  3  km  or  the  energy  ratio  is  lower  than  1.15  no  capture  occurs. 

P0 

Fig  7  19  answers  the  question  about  the  optimal  penetration  speed  of  the 
evader  to  maximize  the  time  of  capture.  The  results  show  that  the  almost 
transonic  speed  is  the  best.  It  is  also  shown  that  for  initial  speed  lower 
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than  250  m/sec  the  evader  cannot  follow  the  optimal  strategy  without  violating 
ground  clearance. 

Fig.  7.20  shows  the  influence  of  thrust  to  weight  ratio  of  the  airplanes. 
It  seems  to  be  a  parameter  of  major  importance  since  the  turning  phase  of  the 
interception  is  relatively  short.  Using  an  example  which  under  nominal 
conditions  has  a  time  of  capture  of  t^  =  216  sec,  it  is  shown  that  a  5% 
improvement  of  pursuer's  thrust  to  weight  ratio  (or  equivalently  5%  degrada¬ 
tion  for  the  evader)  can  reduce  this  capture  time  by  about  30%.  The  slope  of 
the  curve  for  the  nominal  value  indicates  that  navigation  to  the  opposite 
direction  (increasing  thrust  for  the  evader)  does  not  allow  capture. 


TABLE  7.1.  Aircraft  Data 


— 

Aircraft 

A 

B 

W[Kgf] 

12700 

20000 

S[m2l 

27,87 

49.246 

W/S[kgf/m2] 

455.7 

406 

T0/W 

0.86 

0.47 

n 

max 

6 

5 

M 

max 

2 

2.1 

. 

5 


TABLE  7.2.  "A”  AIRCRAFT  THRUST  AND  AERODYNAMIC  DATA 


0.6 

0.8 

1 

1.2 

1.4 

1.6 

1.8 

2 

0 

isa 

11793 

13109 

15060 

15105 

14970 

13562 

12111 

3050 

8618 

9526 

10523 

12247 

12973 

13517 

12292 

11340 

6100 

6260 

7258 

8255 

9571 

10205 

11431 

11340 

10796 

9150 

4173 

5035 

6124 

7258 

8210 

8845 

9117 

mm 

12200 

2540 

3084 

5534 

6033 

6396 

6622 

15250 

1406 

1769 

2268 

2722 

3220 

3493 

3719 

3856 

Fig.  7.1:  Aircraft's  A  thrust  coefficient. 


Aircraft’s  A  parasitic  drag  coefficient. 


Fig. 


Aircraft's  A  maximum  lift  coefficient. 
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Time  history  of  the  bank  angles  (positive  for  left  turn). 


Fig.  7.9:  Time  history  of  azimuth  and  horizontal  line  of  sight  angles. 


Fig.  7.10:  Time  history  of  flight  path  and  vertical  line  of  sight  angles. 
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Effect  of  pursuer’s  initial  alt  it 


18:  Effect  of  pursuer's  initial  specific  energy  variation  on  capture  time. 


19:  Effect  of  evader's  initial  velocity  variation  on  capture 


lip.  7.20.  Kffect  of  thrust  to  weight  ratio  variation  of  both  aircraft  on  th 


SECTION  VIII 


ACCURACY  ASSESSMENT  OF  THE  ZERO-ORDER  FSPT  SOLUTION 


An  important  test  for  the  usefulness  of  any  approximation  is  its  accuracy. 

J  Accuracy  depends  generally  on  satisfaction  of  the  simplifying  assumptions  made 

in  order  to  obtain  an  analytically  treatable  problem.  The  solution  of  the 
medium  range  interception  game  by  the  method  of  forced  singular  perturbations, 

|  presented  m  previous  chapters,  seems  to  demonstrate  that  the  basic  hypothesis 

of  time  scale  separation  of  variables  has  been  a  valid  and  justified  assumption. 
This  time  scale  separation  is  however  not  perfect  and  as  a  consequence  the 
1  accuracy  of  the  zero-order  feedback  approximation  cannot  be  expected  to  be 

absolute. 

The  accuracy  of  any  approximation  has  to  be  established  by  comparison  to 
^  the  results  of  the  exact  solution.  Unfortunately,  there  is  no  operating 

numerical  algorithm  to  solve  the  non-linear  two-point  boundary  layer  problem 
formulated  by  the  necessary  conditions  of  saddle  point  optimality  for  a  3-D 

I 

'  time  optimal  variable  speed  pursuit -evasion  g.imo. 

In  the  last  decade  several  computer  programs  were  developed  in  order  to 
^  satisfy  this  obvious  need  (20-27]  but  most  of  them  are  not  operative.  The  few 

that  are  in  operating  conditions  [Jarmark,  Enjalbar]  are  limited  to  solve  fixed 
time  games.  Moreover,  they  exhibit  serious  difficulties  of  convergence  for 
,  relatively  long  flight  times  as  in  the  case  of  a  medium  range  interception. 

As  a  consequence  of  the  state  of  art  of  reliable  numerical  game  solutions, 
this  report  cannot  include  a  quantitative  assessment  of  accuracy  of  the  zero- 
I  order  FSFT  feedback  solution  for  the  3-D  medium  range  air-to-air  interception 

game.  It  must  be  mentioned,  however,  that  the  one  sided  optimal  control 
version  of  the  same  3-D  medium  range  interception  engagement  (assuming  a  known 

i 
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evader  control  strategy,  as  for  example,  straight  flight)  has  been  also  solved 
by  the  similar  FSPT  approach  [18].  This  optimal  control  approximation  has 
been  numerically  validated  by  comparing  to  the  results  obtained  by  a  very 
accurate  multiple  shooting  algorithm  [19].  For  the  computations,  the  aero¬ 
dynamic  and  propulsion  model  of  aircraft  "B",  in  this  report,  was  used.  The 
results  of  this  comparison  were  presented  in  an  invited  paper  at  the  JACC  1981 
[  1  c, |  The  comparison  showed  an  excellent  qualitative  similarity  in  the  time 
history  of  the  main  variables  and  an  accuracy  of  the  order  of  1 %  (or  better) . 

Very  recently  a  numerical  comparison  was  carried  out  at  NASA  Ames  Research 
Center  between  the  zero-order  FSPT  approximation  and  the  exact  solution  for  a 
variable  speed  horizontal  pursuit-evasion  game.  Results  of  the  comparison  are 
summarized  in  a  forthcoming  document  [28] .  For  this  comparison  the  exact 
solution  was  obtained  by  backward  (retrograde)  integration  of  the  state  and 
adjoint  equations  expressed  in  the  coordinate  frame  of  the  final  line  of  sight. 
This  method  was  proposed  to  obtain  "open-loop  optimal"  extremal  trajectories  in 
several  recent  papers  [29-31].  The  outcome  of  a  systematic  comparison  for  a 
large  set.  of  different  initial  conditions  has  demonstrated  that  for  all  cases 
where  the  geometrical  perturbation  parameter,  defined  by  the  ratio  of  the 
evaders'  minimum  turning  radius  at  the  initial  flight  conditions  to  the  initial 
range  of  separation  is  sufficiently  small  (less  than  0.1)  the  accuracy  of  the 
zero-order  FSPT  solution  is  1%  or  better. 

Although  these  last  results  of  game  comparison  were  obtained  for  a  two- 
dimensional  case  (there  is  no  similar  open  loop  solution  for  the  3-D  game)  they 
are  consistent  with  the  results  of  [19]. 

It  can  be  therefore  asserted  that  the  zero  order  FSPT  approach  is  a  valid 
method  of  analysis  for  the  medium  range  air-to-air  interception  game  (since  in 


this  case  the  conditions  of  geometric  time  scale  separation  are  almost  always 
satisfied)  even  if  presently  no  direct  numerical  assessment  can  be  accomplished. 


The  accuracy  of  1%  or  better,  indicated  in  both  related  comparisons,  seems 
|  to  be  satisfactory  for  all  practical  purposes  such  as  a  parametric  study 

presented  in  Ch.  7  or  for  a  real  time  airborn  implementation.  The  accuracy 
of  regular  state  variable  measurements  and  estimates  of  performance  parameters 
^  does  not  attain  that  level  any  longer.  Consequently,  the  computation  of  first 

and  higher  order  terms  of  the  FSPT  solution  may  not  be  necessary.  Nevertheless, 
since  such  correction  terms  may  be  required  for  some  special  cases,  the 
technique  of  their  derivation,  which  has  been  a  part  of  the  research  program, 
is  presented  in  the  next  chapter. 

! 

» 

I 

I 
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SECTION  IX 


L 

I 

i 

< 


! 

I 
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OFF-LINE  FIRST-ORDFR  CORRECTIONS 


1.  GENERAL  DISCUSSION 

Comparing  the  expressions  for  optimal  control  strategies  of  the  exact  op¬ 
timal  saddle  point  solution  to  the  expressions  of  the  zero-order  multiple  time 
scale  FSPT  approximation, it  can  be  directly  observed  that  the  functional 
structure  (i.e.,  the  dependence  on  the  state  and  the  adjoint  variables)  is 
identical.  The  only  difference  is  that  in  the  FSPT  solution  some  of  the  adjoint 
variables  are  approximated  by  a  simplified  state  feedback  expression, obtained  by 
assuming  a  time  scale  separation  of  the  state  variables.  Improving  the  zero- 
order  feedback  solution  means  to  use  a  better  estimate  for  the  adjoints.  Con¬ 
sequently,  the  control  strategy  will  be  also  improved. 

For  computation  of  the  higher  order  terms  the  method  of  matched  asymptotic 
expansion  can  be  used.  It  has  to  be  however  remembered  that  since  we  deal  with 
a  feedback  solution,  the  state  variables  are  assumed  to  be  accurately  measured. 
Only  the  adjoint  variables  and  the  controls  have  to  be  expanded.  Such  an  ex¬ 
pansion  is  different  from  usual  open  loop  expansions. 

As  it  will  be  shown  in  this  chapter,  the  computation  of  the  first  order 
correction  terms  has  to  be  carried  out  by  integration  along  the  previously  ob¬ 
tained  zero-order  trajectory.  It  is  clearly  an  "off-line"  operation  and  cannot 
be  considered  as  a  part  of  a  real  time  implementation.  It  has  however  its 
engineering  value  if  off-line  computation  is  not  prohibitive  and  if  additional 
accuracy  is  required. 

In  order  to  demonstrate  the  complexity  involved  in  such  operation,  in  the 
following  parts  of  this  chapter  the  example  of  horizontal  medium  range  inter¬ 
ception  (solved  to  zero-order  in  [3])  is  selected.  In  the  sequel  the  non- 
dimensional  notation  of  [3l  will  be  used. 
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LIST  OF  SYMBOLS  FOR  THE  NON-DIMENSIONAL  HORIZONTAL  GAME 


-  Hamiltonian 

-  Components  of  the  adjoint  vector 


-  Mach  dependent  coefficient 

-  Aircraft  turning  radius 

-  Separation  distance 

-  Normalized  separation  distance 

-  Time 


min 


-  Normalized  time 


min 


-  Aircraft  true  airspeed 


Normalized  velocity 


-  Ratio  of  minimum  turning  radii 


P 

(v  A  min\ 

v  a*  ) 


E  . 
min 


-  Small  parameter 

-  Normalized  turning  rate 


-  Line  of  sight  angle 

-  Normalized  acceleration  (A  =  $  -  Nn^) 


-  Stretched  time  in  the  inner  boundary  layer  (a  =  t /e^) 

-  Stretched  time  in  the  outer  boundary  layer  (x  =  t/e) 

-  Aircraft  straight  line  normalized  acceleration 

-  Velocity  vector  angle 

(  A  VpminN 

-  Ratio  of  corner  velocities  I  Q  = - _) 


SUBSCRIPTS 


C  -  Comer  conditions 

E  -  Evader 

P  -  Pursuer 

h  -  Horizontal  plane  value 

0  -  Zero -order  value 

1  -  First-order  value 

b  -  Value  at  termination  of  the  boundary  layer 

SUPERSCRIPTS 

o  -  Reduced  game 

i  -  Outer  boundary  layer 

j  -  Inner  boundary  layer 

c  -  Composite  variable 
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3. 


FIRST  ORDER  EQUATIONS 


~  =  -vp  cosC(pp  -8h)  +  fiv£  sin((pE  -  8h) 

-Jl  .  I  [-VP  sin(«>p  -eh)  *  nvE  sin«PE  -9h)] 
dt 


dv 

e  — r  -  <f>P  -  Np^np  +  enp  ) 
dt  r  0  1 


dvp  r  2 " 

e  -4  =  rn  U  -  n  (nE  +  enE  ) 

dt  •  t  to  fci 


d<P 

_  =  n  +  en 
dt  0  1 


2  d<J)E 

£  _£  =  TQ 

dt 


n  +  en 
h0  1J 


The  Hamiltonian  is: 


(j  .CJ  )^.(j  .  CJ  )tl 

V  Vn  vD  )  A~  \  vc  vc  )  A~ 


\  V  *  *  %'  * 


*  (\  *  e\X\  *  c\)  *  (j<pe„  *  s,)  rnK  *  cn0 


E  *  enE 

p  -T'p  /  \ 

c0  C1  u 


Since  time  does  not  appear  explicitly  in  the  equations  of  motion: 


.If  =  3fQ  ♦  e3f  =  0 


where : 


(9.1) 

(9.2) 

(9.3) 

(9.4) 

(9.5) 

(9.6) 


(9.7) 


(9.8) 
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(9.9) 


3f  = 

,  T  dr 

1  +  J  - 

0 

r0  dt 

J  | 

[^)  w 
^  dt'  >0  vl 

VE 

&) 
v  dt  7< 

•(-) 

D  Xdt70 

rdvp\ 

(,dvp\  , 

^dVp\ 

*  (lT'o 

dv 


,  de  .  / 
%  *  \  ( 


n  +  j  ra  n„  =o 
P  0  po  %  Eo 


/  ve\  /uvp\ 

I - 1  ,  ( - 1  —  zero-order  terms  of  the  aircraft  accelerations 

vdt70  vdt70 


t dv . 


('dv. 


dt 


+  n„  *  J,„  n,  +  TO  .1  n  +  J  ri  =  0 


r 


p 

po  1 


ip  1 

I’j  0 


Lo  1 


1 


Jtp  E  I 

E}  0J 


The  adjoint  equations  are: 


dJ  dJ 

ro 

■  ■  —  +  c 
dt  dt 


r-  =  (v  *  E-0  4 1  -vp  sin(<pj>  -v  *  ""e  sin(,,E  -  v 

t  u  l  r 


"ft  dJ0 

- -  +  e  — 

dt  dt 


(Jr0  *  EJri)  [VP  sint<pp  -  8h)  -  [JvE  sin(«)E-8hl 


*  (\  *  eJ(ii) ;  ['Vp  cos(,®p  - eh>  *  ^E  cosW>E  -  v] 
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(9.10) 


(9.11) 


(9.12) 


(9.13) 


(9.14) 


dj  dJ 

Vpo 

e  _ _  ♦  e  - : 


— -  =  ( J  *  ej  )  cos  (ipn  - 

t  V  r0  V  P 


V  + 


*  (Je„  *  "ej  |  sintipp  ■  V  -  (Jv.  4  ^v.  ) 


_  /  °V_ 

Pj'  P 


(9.15) 


dj  dJ 

VE  VE 

£  -Jo  +  £  _Jl  .  -Q 


Ks  *  s) 


cos(<P_  -  6.  )  ♦ 
t  n 


+  (jn  ♦  £Jn  )  -  sin(<P  -  0.  )  -  (j  +  eJ  )  TQ 

V  e0  Vi  E  h  V  vE  vE  y  ^  J 


(9.16) 


dJiPp  ^^(Pp 

e2  - 5.  +  £  - L  -  v  +  £j  )  sin(y>  -0  )  + 

di  di  p  L  v  ro  v  p  h 


4  (Je  4  m0i) 

0  It 


(9.17) 


dd<p  dJcp 

2  E0  E 

E.  ■■■  ■  +  £  1  ■  ■ 


— “  =  [fj  +  eJ  1  sinOP-  -0.)  - 

t  E  lV  r0  V  E  h 


[Je0  4  £,e1)r“»<V<lh 


(9.18) 


4.  THE  REDUCED  GAME 


During  the  reduced  order  game: 


(9.19) 


8A 


MICROCOPY  RESOLUTION  TEST  CHART 

national  buacau  of  standards  —  IVCS  ~  A 


From  (9.6),  (9.5) 


El 


1?  -  "e  ■  0 


Substituting  (9.20),  (9.19)  into  (9.12) 


EH? 


o  ,°  / dr\° 

1 '  \  W 


& 


However,  since  { — j  t  0 
Mt' 


J°  =  0 


(9.20) 


(9.21) 


(9.22) 


Optimality  condition  of  the  reduced  game  (Eq.  (69)  of  [3])  is: 


S 


v 

vp  sin(ipD  -0J  -  —  cos (<PD  -  0V )  =  0 


rP  h  0 


P  h 


(9.23) 


i 


During  the  reduced  game: 


< 


=  0 


(9.24) 


Therefore: 


J°  =  J°  =0 

eo  ei 


(9.25) 


Finally,  from  (9.9),  (9.19),  we  obtain: 


0  vj  -  Qv° 


(9.26) 
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5.  THE  OUTER  BOUNDARY  LAYER 

Using  the  stretching  transformation  T  =  t/e,  the  following  adjoint  equations  are 
obtained: 


dJ 


dJ 


+  £ 


dT 


(JeQ  *  CJej)  4  [-vp  -  V  *  f!''Sin«PE-8h)] 


(9.27) 


dJ 


dx 


> e  -ar  ■ 6  (Ji0  *  eJij)  [VP sin'^  •  V  - 

■a  vE  sin(<PE  -eh)]  *  £  (jJ  *  jj  )  7  [-vp  cos(4.  -e  ) 
j  '0  V  r  L 


+  cos(^E  -  0h) 


(9.28) 


dJ 


dJ1 


dT 


c  -ar1  ■  (Ji0  *  "ij)  cosWp  -v  * 
(Je0*«e1)r-"^-V  -  (4p  *  <p  )  % 


(9 . 29) 


dJ 


VE, 


dJ1 


dT 


+  £  - 


dT 


-a  [(Jr0  *  eJij) 


(4  *  r  -  v]  -  (JvE  *  «t„  )  rn  5^ 


<p  <P 

Po  ,  e 

"Hi  “3t 


vD  [-(j1  +  eJ1  ) 

p  L  V  rQ  xx) 

(Je0  *  £,S1):“!l*p-v] 


sin(<Pp  -6h)  + 


(9.30) 


(9.31) 


i 
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e  ~3t  +  £ 


+  £  itr1}  -  ([Ji0  +  ej1t)  siT,(<PE  -V 

■  (J0O  *  “Sj)  j  cosCl°E  '  v] 


From  the  zero-order  terms  of  (9.27)  we  have: 


__2.  =  0  «  J*  =  const, 
dx  rQ 


(9.32) 


(9.33) 


ri  Ti  If 

it  =  J0  ~2  [V 
Or  L 


p  sin(<Pp  -eh)  +  sin(q^  -0h)j 


and  from  the  zero-order  terms  of  (9.28): 


(9.34) 


=  0  «-»  J  =  const . 
T  90 


(9.35) 


0  r  -i 

IT*  JrQ  l*vP  ‘‘"‘W  *  nvE  slnlWJ  * 

*  J60  f  ["Vp  C°St<5'’  'V  *  ^  C0SWE  -  V] 


(9.36) 


Using  (9.33)  and  the  matching  condition: 


lim  J1  =  J°  =  - 

x-nb  r0  r0  t-0  v°  - 


(9.37) 


From  (9.35)  and  the  matching  condition: 


J%  -  ’  ° 


(9.38) 


From  (9.34)  and  (9.38): 


I 
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=  const. 


(9.39) 


This  integral  can  be  computed  only  if  the  zero-order  trajectory  (i.e.,  vp(T),  4>p(T), 
0h(T)»  ve(T)»  Ve(t))  are  known. 

At  this  point  it  is  necessary  to  show  that  in  this  boundary  layer: 
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Ap  -  *p  +  °Ce  ) 


(9.46) 


From  (9.5) : 


d*P  i  i 

£  ar =  %  +  enPl 


(9.47) 


It  is  clear  that: 


n;  =  0 

0 


(9.48) 


Substituting  (9.46)  into: 


ap  ’  *P  -  np  ("p0  *  "pj 


(9.49) 


we  get 


(9.50) 


First-order  terms  of  (9.29)  are: 


P1  i  1  i  /8V 

\  Tsi”<VV  -  (*y 


(9.51) 


(9.51)  can  be  integrated,  now,  starting  from  the  terminal  surface  and  using  the 
boundary  condition: 


Ix  Tv 

v„  v„ 

p  p  ■ 

1  0 


(9.52) 


This  integration  can  be  done  only  if  the  zero-order  trajectory  is  known. 
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~6.  THE  INNER  BOUNDARY  LAYER 

2 

Using  the  stretching  transformation  a  =  t/e  ,  the  following  adjoint  equations  are 
obtained: 


(9.53) 


<“  i  dJe 

0  +  -  1 

do  do 


«•“!,)’  (»)  w.  •"!,)  KS 


(9.54) 


From  (9.53),  (9.54),  the  following  system  of  equations  is  obtained: 


0  ~  =  J1 


ro  ro  T=0  ro 


(9.55) 


,o  r. 


(9.56) 


-sr-0--*. 


(9.57) 


Tb 

-  0  -  jj  =  jJj  0  3  J°  J  [”VP  sin(Vp-eh)  +«vE  sin  (<PE-eh)  ]dr  (9.58) 


adjoint  equation  is: 
VP 


!o  +  e  ~3o 


To"1  -  E  (Jr0  *  Ui,)  C0S(»P  '  V  * 

-£(Jip  *Mip) 


(9.59) 


The  first  order  equation  (using  (9.55)  and  (9.57))  is: 
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dJJ 

VP 

So1-  Jr0  C0S<W 


-  J 


(UP 


\  H 


•l 


(9.60) 


2  2 
Ap  ■  <frp  -  Np  (nj  +  2enJ  nj  +  e2nj  ) 
k  r  r  \  *o  *0  *1  V 


(9.61) 


when  G  =  const. 


©0  ■  © 


(9.62) 


Integrating  (9.60)  along  the  zero-order  trajectory: 


Jv  =J  cosC<PP  '  V  d° 
?1  0J  0 


•  f->°  (?)* 

o  VP0  V3V 


Boundary  condition: 


(9.63) 


0=0, 


T=0 


°b  *  f°b  O  /3*p\ 

cos(tpn-6.)  do  -  J°  (w — )  do 

T.o  0'  -0  p  h  VP0  V V 


f  b 

=  J  J°, 


(9.64) 


(9.65) 


=  J1 

VP  VP 
*1  1 


|“j°  cos(Wp-9h)d0.  j“j°  (©*>  (9.66) 


T°0 


0J  0 


0J  'P  'W'P' 


J?.  +  J1.  -  J1. 


T=0 


(9.67) 


The  inner  boundary  layer  Hamiltonian: 

**  .  i  .jj  ;  ♦  jj (  «  .  Aj  nj  ♦  jJ  rnA^  *  4  ran* 

P  P  E  E 


(9.68) 
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Let  us  calculate  the  first  order  terms  of  the  following  expression: 


i  •  c « > 


2 

(4  4)  *  4  (*p  -  vjp  )  -  ^p  4  4  4 

\  Vp  vp  \  V  *0'  r  Vp  r0  ri 


The  optimality  condition  is: 


2np  (4  *  <  X4#  *  *4.)  ■  <  *  £4P 


The  first  order  terms  satisfy: 


2Np  (4  4  *  4  40)  ■  4P 


Multiplying  both  sides  by  we  obtain: 

0 


-2np  4  4  4  ■  -4  4  *  2np  4  4 

p0  1  °  pi  °  pi  ° 


Substituting  (9.73)  into  (9.71) 


Ji  "i  * 

vPj  0 


*  2NP  Jv  np  *  Jv  \*T  *  Vp  )  -  Ji 

VP  0  VP  v  *V  VP  0 


Substituting  (9.75)  into  (9.69)  we  get: 


*i  =  Je  §  *  JJ.  (*p  *  np<)  *  -'i.  np,  * 
1  U  0 

.  ra  jj  (*  .  N  p|  )  ♦  rn  =  o 

E,'  O'  1 


The  first  order  correction  term  for  the  pursuer's  turn  rate  control  is 


.  I A  Pl^JVV  1 ,  (*  .  -  -tf  jJ.  \ 


:  2\  , 


Similar  expression  can  be  derived  for  the  evader. 

Finally  it  can  be  observed  that  the  entire  right-hand  side  of  (9.77)  is  derived 
from  the  zero-order  trajectory. 

Jo  is  calculated  from  (9.58)  and  J"*  from  (9.64)  and  (9.51). 

1  Vpi 

Both  expressions  require  integration  along  the  zero-order  solution  which  is  an 
off-line  operation. 


7 .  CONCLUSIONS 


Due  to  the  inherent  complexity  of  the  approach  outlined  in  this  chapter,  even  for 
the  relatively  simple  2-D  (horizontal)  case,  it  seems  that  other  methods  for  approving 
the  zero-order  FSPT  approximation  need  to  be  investigated. 

An  example  for  such  an  improvement  technique,  motivated  by  a  recent  work  [32]  is 
presented  in  the  next  chapter. 
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SECTION  X 


IMPROVEMENT  OF  THE  ZERO-ORDER  GAME  SOLUTION  BY  PARAMETER  OPTIMIZATION 

1.  GENERAL  DISCUSSION 

As  mentioned  in  Chapter  9,th*  functional  structure  of  the  suboptimal  control 
strategies  (i.e.,  the  dependence  on  the  state  and  the  adjoint  variables)  is 
identical  to  the  optimal  control  strategies. 

The  difference  is  that  the  adjoint  variables  being  based  on  the  time  scale 
separation  assumption,  are  approximate. 

One  way  to  improve  the  approximation  of  zero-order  is  to  compute  higher  order 
terms  as  demonstrated  in  the  previous  chapter.  However,  the  complexity  of  this 
approach  renders  it  a  hardly  practical  engineering  tool . 

In  this  chapter  a  different  approach  is  presented.  Motivated  by  a  recently 
reported  investigation  by  ONERA  [32]  on  approximating  a  game  saddle-point  solution 
by  results  of  respective  minmax  and  maxmin  optimization  it  is  proposed  to  optimize 
parameters  in  the  feedback  expressions  of  the  FSPT  control  strategies. 

The  most  important  parameters  in  this  respect  are  the  maximum  velocities  of  the 
respective  airplanes,  determined  for  each  of  them  by 


V 

max 


=  arg  max 
E,h 


[T, 


max 


-  D„ 


*Di> 


(10.1) 


In  the  reduced  order  game  it  has  been  assumed  that  both  aircraft  reach  this 
maximum  velocity  before  termination,  i.e., 

VE  =  W  =  VE  .  v?  -  vp  -  W  =  Vp  (10.2) 

max  max 

In  reality  these  conditions  are  rarely  satisfied.  Consequently  the  value  of 

-  l/cv“  -  V°)  .  1/Ky^  -  (vE)max)  (10.3) 

is  not  correct,  violating  the  transversality  condition 


5f(tf)  =  0 


(10.4) 
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Using  the  uncorrect  values  for  Vp(tf)  and  Vg(tf)  leads  to  further  violations 
of  other  transversal ity  conditions  as 

AE  =  XE  =  °  (10.5) 

P  E 

resulting  in  a  deviation  from  optimality. 

An  educated  guess  of  the  actual  final  velocities  of  the  opponent  airplanes  (which 
are  not  known  a  priori)  Vp(tf)  and  V£(tf),  and  their  use  in  the  zero-order  feed¬ 
back  control  expressions  instead  of  t^p^max  and  ^p-^max’  s^ou^d  provide  an  improved 
approximation  for  the  game  solution.  Such  a  correction  should  also  lead  to  improve 
the  performance  of  the  individual  aircraft. 

For  purposes  of  performance  assessment  by  off-line  computation  this  correction 
can  be  easily  carried  out  by  successive  approximation  starting  with  the  zero-order 
feedback  solution.  At  each  step  the  values  of  V°  and  V°  are  approximated  by  the 
actual  final  values  obtained  from  the  previous  step.  The  convergence  of  this  process 
is  very  fast.  Experience  gained  using  the  examples  of  Chapter  7  indicates  that  only 
2  or  3  additional  integrations  are  necessary  and  the  improvement  can  be  of  importance. 

For  a  real  time  (on-line)  airborne  implementation  (one  of  the  most  attractive 
features  of  the  FSPT  approach),  simple  algorithms  can  be  developed,  based  on  range 
measurements  and  estimations  of  own  and  opponent  acceleration  capabilities,  to  predict 
approximative ly  V(t^)  of  each  airplane.  The  type  of  feasible  algorithm  depends  on 
the  available  measurements  onboard  the  airplane.  Further  exploration  and  assessment 
of  such  algorithms  are  out  of  the  scope  of  the  reported  research  activity. 


'2.  THE  PARAMETER  OPTIMIZATION  PROCESS 

In  this  section  the  effect  of  varying  the  value  of  V°  and  V°  on  the  outcome  of 
the  game  is  demonstrated.  In  order  to  analyse  the  respective  performance  improvements 
and  verify  that  optimal  value  of  V°  the  approach  of  minmax  (maxmin)  parameter  optimi¬ 
zation,  motivated  by  [32],  was  selected.  The  example  for  this  analysis  was  the  one 
used  in  Chapter  7  with  the  following  set  of  initial  conditions: 

R  =  20.000  m,  ip  =  0°,-Xd  «  0°,  xc  =  180°,  V  =  280  m/sec,  V  =  370  m/sec, 

u  o  PQ  Eq  Pq 

hD  =  4000  m,  hp  =  2000  m,  and  with  R_  =  4000  m. 

*0  u  h 

The  results  of  the  uncorrected  zero-order  solution  are  summarized  in  Table  10.1. 
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The  large  discrepancies  (especially  for  the  evader)  indicate  that  this  particular 
zero-order  game  solution  may  not  be  a  good  approximation  of  an  optimal  game. 

To  evaluate  the  influence  of  V°  on  the  individual  aircraft  performance  first  the 
evader's  V°  was  varied  as  a  parameter,  keeping  Vp  =  (Vp)max  =  555  m/sec  constant. 

Results  of  these  computations  are  depicted  by  a  solid  line  ofi  Fig.  10.1  as  a  function 
of  the  speed  ratio  (v°/vmax) •  The  same  excercise  was  repeated  by  varying  Vp  and 

keeping  V°  =  (V  )  =  604  m/sec.  This  results  are  shown  in  Fig.  10.1  by  a  dashed 

fc  b  max 

line.  In  Fig.  10.2  the  ratio  of  the  actual  final  velocities  (V(t-)/(V  ))  are 

6  o  max" 

plotted  as  a  function  of  (V  /V  )  for  each  case. 

max 

These  results  confirm  the  expected  performance  improvements  as  V°  -*•  V(tf), 

amounting  to  an  individual  sensitivity,  defined  as  the  ratio  of  the  relative  change 

in  capture  time  to  the  relative  change  in  the  value  of  own  V°  (while  the  opponent's 

V°  =  V  ) . 
max 


i_W_ 

'  K'W’  V?  -  (V„ax 


=  -0.33 


(10.6) 


(Atf  /tf) 

K/W]  VE  =  (Vmax 


=  0.13 


(10.7) 


As  expected  by  optimizing  the  parameter  V°  for  the  evader  (the  maximizer)  the 
capture  time  increases  and  conversely  by  varying  Vp  to  the  right  direction  the 
capture  time  is  decreasing. 

The  result  of  a  simultaneous  parameter  optimization  is  not  obvious  from  the 
beginning,  but  due  to  the  higher  sensitivity  of  the  evader  it  can  be  expected  that 
the  corrected  game  solution  would  favourably  shift  to  its  direction. 
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Indeed  the  outcome  of  the  engagement  with  converged  values  of  Vp  =  Vp(t^.)  =  507  m/sec 
and  V„  =  V„(t,.)  =  476  m/sec  resulted  in  t^  =  163  sec  an  improvement  of  8.5%  for  the 

C  t  X  “X 

evader  compared  to  the  uncorrected  zero-order  estimate  of  149  sec. 


3.  A  QUALITATIVE  ANALYSIS 

In  Fig.  10.3  and  10.4  the  velocity  profiles  of  the  evader  and  the  pursuer  for  the 
parameter  optimized  engagement  are  respectively  compared  to  their  profiles  in  the 
uncorrected  game. 

Since  the  initial  conditions  of  the  engagement  are  very  unfavourable  for  the 
evader  (y  =  180°)  its  turning  performance  at  the  initial  phase  is  predominant.  The 

0  o 

corrected  lower  value  of  Vp  (see  Eq.  5.103  and  6.12  )  requires  higher  turning  rate 
and  consequently  a  smaller  gain,  or  even  some  loss,  of  specific  energy.  Inspite  of 
the  fact  that  the  game  terminates  later,  the  final  speed  of  the  evader  is  less  than  in 
the  uncorrected  version.  The  relative  advantage  gained  by  this  faster  initial  turn, 
and  demonstrated  by  the  increased  "capture  time",  is  illustrated  in  Fig.  10.5  where 

the  histories  of  the  relative  range  of  separation  of  the  engagements  are  compared. 

The  variations  in  the  speed  profile  of  the  pursuer  in  Fig.  10.4  are  much  less 
impressive  for  the  following  reasons:  (i)  the  pursuer  starts  with  the  most  favourable 
initial  geometry  (y  =  40  and  almost  does  not  have  to  turn  in  the  horizontal  plane; 

P°  o°  o 

(ii)  the  difference  AVp  given  in  Table  10.1  is  much  smaller  than  AVg.  Consequently 

due  to  the  longer  flight  time  the  final  velocity  is  higher  than  in  the  uncorrected 

game. 

The  comparison  also  indicates  that  the  relatively  large  improvement  is  very  probably 
due  to  the  asymmetry  of  initial  conditions  of  the  example  selected  and  to  the  fact  that 

the  evader  has  a  higher  maximum  velocity  than  the  pursuer,  i.e.,  that  reduced  order  game 

does  not  predict  capture. 

The  additional  insight  gained  from  this  comparison  and  the  demonstrated  improvement, 
indicate  the  potential  of  this  simple  technique  and  the  need  for  further,  more  detailed 
investigation  in  this  direction. 
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variation  on  the  capture  time. 


max 


final  velocities. 


Fig.  10.3:  Uvader's  velocity  profile  for  zero-order  and  corrected  games. 


10.4:  Pursuer's  velocity  profile  for  zero-order  and  corrected  games. 


R  l  km  ] 


histories  of  separation  distance  in  zeTO-ordeT  and  corrected  Rames . 


SECTION  XI 


CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

From  the  results  presented  in  this  report  it  can  be  first  of  all  be 
concluded  that  modern  air  combat  can  be  meaningfully  investigated  by  an 
analytical  method.  An  air  combat  problem  of  considerable  importance:  the 
three-dimensional  medium  range  air  to  air  interception,  with  advanced 
guided  missiles,  was  successfully  analyzed  using  realistic  airplane  models, 
by  combining  the  technique  of  forced  singular  perturbation  and  differential 
game  theory.  This  methodology  has  provided  a  zero-order  feedback  control 
strategy,  closely  approximating  the  exact  solution.  The  attractiveness  of 
the  feedback  solution  for  airborne  implementation  cannot  be  overemphasized. 
(This  aspect  of  the  solution  merits  some  additional  effort,  outlined  in  the 
sequel).  The  feedback  solution  can  serve  also  as  a  very  efficient  tool  for 
rapid  but  systematical  assessment  of  several  operational  problems  of 
interest  as  demonstrated  by  the  analysis  of  an  imaginary  air  defence 
scenario. 

The  numerical  accuracy  of  the  zero-order  FSPT  solution  for  a 
three-dimensional  interception  game  could  not  be  directly  assessed  due  to 
the  non-existence  of  any  operating  numerical  algorithm  to  provide  the  exact 
game  solution  for  comparison.  There  is  however  a  very  strong 
"circumstantial  evidence",  based  on  similar  but  not  identical  comparisons, 
that  the  zero-order  FSPT  solution  has  a  sufficient  level  of  accuracy  for 
medium  range  air  combat  analysis.  Consequently  higher  order  correction 
terms  in  many  cases  may  not  be  required. 
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Nevertheless,  the  methodology  to  obtain  the  first  (and  higher)  order  terms 
of  the  l:SPT  expansion  has  been  developed,  indicating  that  such  correction  can  be 
carried  out  only  by  an  off-line  computation.  The  considerable  computational 
complexity,  demonstrated  by  the  relatively  simple  horizontal  example,  makes  the 
usefulness  of  such  correction  questionable.  As  an  alternative,  a  technique  of 
parameter  optimization  is  proposed  and  its  potential  for  improvements  is 
demonstrated.  On-board  implementation  of  this  technique  seems  possible  but 
requires,  however,  further  investigation. 

Another  problem  addressed  in  the  research  program  was  the  solution  of 
"terminal  boundary  layer".  It  was  found  that  terminal  boundary  layer  solutions, 
though  they  may  exist  and  computed  in  open  loop,  cannot  be  implemented  in  a 
feedback  form.  Moreover,  it  has  been  demonstrated  that  solution  of  terminal 
boundary  layers  can  be  avoided  by  using  a  judiciously  formulated  dual  model¬ 
ling  (see  Chapter  5) .  This  approach  enables  to  keep  the  zero-order  control 
strategies  in  a  pure  feedback  form  and  make  airborne  implementation  of  these 
strategies  a  simple  task. 

The  topic  of  airborne  implementation,  being  the  ultimate  objective  of  an 
applicable  analysis,  merits  further  elaboration. 

2 .  AIRBORNE  IMPLEMENTATION 

As  mentioned  earlier,  the  feedback  form  of  the  FSPT  solution  makes  it  an 
extremely  attractive  candidate  for  implementation  onboard  future  advanced 
aircraft.  There  are  several  ongoing  experimental  programs  exploring  concepts 
and  technologies  for  advanced  integrated  fire  and  flight  control  systems. 
Incorporating  elements  of  optimal  air  combat  strategy  in  such  a  system  seems 
to  be  an  obvious  step.  The  simplicity  of  the  FSPT  control  algorithm  and  its 
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accuracy  for  medium  range  air  combat  lead  to  strongly  recommend  its  airborne 
implementation. 

In  order  to  provide  a  preliminary  tool  for  evaluating  such  effort,  the 
requirements  for  the  airborne  implementation  are  briefly  obtained.  It  has  to 
be  noted  that  this  outline  is  by  itself  an  introductory  one  and  does  not 
pretend  to  substitute  a  more  profound  assessment.  For  airborne  implementation 
the  following  requirements  can  be  listed: 

a)  storage  of  aircraft  performance  data  and  results  of  some  off-line  computa¬ 
tions; 

b)  measurement,  filtering  and  data  reduction  of  the  state  variables  involved; 

c)  computation  and  synthesis  of  the  control  algorithm; 

d)  display  of  the  required  control  strategy  or  its  application  as  an  input 
of  the  flight  control  system. 

2.1  Storage  Requirement 

The  largest  storage  requirement  relates  to  the  aircraft  aerodynamic  and 
propulsive  model  as  a  function  of  speed,,  altitude,  and  plane  configuration. 
Either  tabulated  data,  requiring  some  interpolation  scheme  o'"  a  polynomial 
approximation  can  be  used,  as  it  has  been  done  for  a  ground  base  simulation. 
Such  data  can  be  also  used,  and  in  some  airplanes  it  has  been  already  done, 
for  other  purposes  such  as  automatic  weapon  delivery,  flight  path  managements, 
etc.  The  additional  data  to  be  stored,  useful  mostly  for  medium  range  inter¬ 
ception,  is  the  precomputed  "optimal"  altitude  for  best  acceleration  "hj",  as 
a  function  of  specific  energy.  The  eventual  constraints  of  the  flight 
envelop  are  included  in  any  modem  flight  control  system. 
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2.2  Measurement  Requirements 

In  order  to  compute  the  FSPT  control  strategy  the  knowledge  of  the  following 
variables  is  needed: 

a.  velocity  vector,  magnitude  and  direction  in  a  fixed  (inertial)  coordinate 
system, 

b.  altitude  (pressure  altitude  is  sufficient), 

c.  angle  of  attack, 

d .  bank  ang 1 e , 

e.  radar  range  to  the  target, 

f.  line  of  sight  orientation  in  aircraft  body  axes, 

g.  aircraft  actual  weight, 

h .  weapon  system  envelope . 

It  is  easy  to  see,  without  further  detail,  that  the  required  measurements  do 
exist  in  every  modem  military  airplane.  If  line  of  sight  angles  are  beyond  the 
usual  range  of  airborne  radar  gimbals,  this  information  has  to  be  either  a  manual 
input  of  the  pilot  or  a  data  linked  message  by  a  ground  based  controller  (giving 
the  interception  command) . 

2.3  Computations 

Filtering,  transformation  of  coordinates  and  other  phases  of  d-ta  reduction 
may  require  some  computational  load  as  for  any  other  automated  system.  It  is 
strongly  dependent  on  the  type  of  the  airplane  considered  for  implementation. 

The  computation  of  the  control  algorithm  is  extremely  simple.  "Hiree  different 
types  of  optimal  load  factors  (n^,nv,ny)  have  to  be  computed.  The  number  of 
arithmetical  operations  in  each  computation  is  relatively  small  depending 


slightly  on  the  type  of  stored  data.  Based  on  these  three  values  the  synthesis 
of  the  required  optimal  load  factor  and  bank  angle  is  obtained. 

2.4  Display 

The  outputs  of  the  control  algorithm  are  the  required  values  of  variables 
already  displayed  in  almost  every  aircraft.  It  is  very  simple  to  use  the 
already  existing  symbols  of  a  HUD  or  other  type  of  display,  by  selecting  a 
proper  mode  of  operation  (i.e,,  medium  range  air-to-air)  in  order  to  indicate 
to  the  pilot  the  required  control  action  minimizing  the  time  for  air-to-air 
weapon  deployment.  Simultaneously,  the  signals  of  the  required  maneuver  can 
!e  used  as  inputs  to  an  automatic  flight  control  system. 


3 .  RECOMMENDATIONS 

Ill  spite  of  the  preliminary  nature  of  the  airborne  implementation  scheme  it 
is  clear  that  it  does  not  present  any  paiticular  difficulty.  Moreover,  very 
probably  it  can  be  carried  out  without  substantial  complication  and  additional 
cost  in  a  planned  integrated  fire  and  flight  control  system.  Therefore  it  is 
strongly  recommended  to  make  the  necessary  steps  that  the  airborne  implementation 
of  the  FSPT  algorithm  for  medium  range  air  combat  would  be  included  in  a  future 
experimental  flight  test  program.  The  real  merits  of  the  algorithm  can  be  fully 
evaluated  based  on  the  results  of  flight  experience. 

further  investigation  of  simple  correction  techniques  of  the  basic  zero-order 
FSPT  solution,  which  seem  to  lead  to  important  improvements,  are  also  recommended. 
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